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Power  requirements  for  future  space  satellites  are  expected  to  increase  significantly 
as  neWN capab ilit ies  for  communications,  data  processing  and  surveillance  are  imple¬ 
mented  ^“More  effective  and  efficient  thermal  transport  techniques  will  be  needed  for 
heat  rejection  from  equipment  on  satellites.  Circulating  two-phase  fluid  loops 

have  been  suggested  and  laboratory  tested  for  possible  application  in  the  above  areas. 

In  comparison  to  a  single-phase  loop,  the  two-phase  system  operates  at  considerably 
smaller  flow  rates  and  maintains  a  tighter  temperature  control  with  higher  heat  transfer 
coefficients.  However,  the  two-phase  fluid  flow  regimes,  pressure  gradients  and  heat 
transfer  coefficients  must  be  evaluated  for  application  in  the  weightless  environment 
of  an  orbiting  satellite. 

This  projectTwas  funded  by  the  Air  Force  Wright  Aeronautical  Laboratories  as  a  Phase  II 
SBTR  program  to^stud-v  two-phase  flow  behavior  under  zero  gravity  conditions.  The  overall 
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objectives  of  this  study  were  to  generate  a  data  base  for  two-phase  pressure  drop 
and  the  void-quality  relationship  under  simulated  zero  gravity  conditions  and  to 
develop  analytical  models  to  predict  these  parameters  for  bubbly  and  annular  flow. 


The  simulation  of  zero  gravity  two-phase  flow  was  achieved  by  using  two  immiscible 
liquids  with  equal  densities  to  eliminate  the  buoyancy  component.  Although  this 
approach  does  not  eliminate  the  gravity  effects,  it  provides  a  representation  for  void 
distribution  in  the  absence  of  gravity..  Water  was  used  as  one  of  the  working  fluids 
and  simulated  the  gaseous  phase,  whiLe'che  second  fluid  was  selected  to  simulate  the 
liquid  phase.  The  property  ratios  of  the  two  liquids  were  in  the  range  of  liquid/gas 
flows.  The  tests  were  performed  under  steady  state  fully  developed  turbulent  flow 
conditions,  since  it  is  believed  that  for  most  practical  applications  the  two-phase 
heat  transfer  loops  operate  within  the  turbulent  flow  range.  The  major  parameters  are 
pressure  drop  and  cross  sectional  average  void  fraction,  and  the  test  variables  con¬ 
sist  of  total  flow  rate,  quality,  and  fluid  properties. 


>The  modeling  effort  is  limited  to  developing  relations  for  the  two-phase  friction 
multiplier  and  void-quality  relation  under  bubbly  and  annular  flow  conditions.  The 
bubbly  flow  model  is  based  on  the  assumption  of  local  homogeneous  conditions  between 
the  phases  but  allows  for  void  distribution  in  the  radial  direction.  Separated  flow 
conservation  equations  are  used,  and  single-phase  turbulent  flow  eddy  diffusivity 
relations  are  employed,  yin  addition,  a  single  particle  lift  force  model  is  utilized 
to  accounc  for  radial  forces  on''  the  bubbles  within  the  shear  layer.  A  tentative  model 
is  proposed  for  bubbly  flow  with  nucleation  at  the  wall.  This  aproach  uses  the  model 
developed  for  adiabatic  flow,  but  the  boundary  condition  is  modified  to  account  for 
the  presence  of  voids  at  the  wall.  The  annular  flow  model  is  based  on  the  triangular 
approach  of  Hewitt.  The  interfacial  shear  and  entrainment  relations  are  empirically 
determined  from  the  results  of  the  experimental  effort  with  equal  density  liquids. 


It  should  be  noted  that  in  the  present  analytical  effort  no  attempt  was  made  to  model 
the  heat  transfer  characteristics  of  two-phase  flow  in  zero  gravity.  It  is  believed 
that  the  mechanism  of  heat  transfer  will  not  be  affected  by  the  absence  of  gravity. 
However,  the  heat  transfer  rates  are  influenced  by  void  distribution  within  the  two- 
phase  flow.  The  earth  gravity  models  are  generally  based  on  flow  quality  and  therefore 
cannot  be  applied  to  reduced  gravity  conditions.  If  one  can  incorporate  the  differences 
in  void-quality  relation  into  the  earth  gravity  correlations,  they  can  be  used  under 
reduced  gravity  conditions. 
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Section  1 


INTRODUCTION 


The  power  requirements  for  future  spacecraft  and  satellites  are  expected 
to  increase  as  they  grow  in  sophistication  and  capability.  Due  to  the 
new  communications,  data  processing,  and  surveillance  techniques,  the 
near  future  satellites  are  planned  to  operate  at  power  levels  of  five  to 
ten  kilowatts.  With  the  future  revolutionary  military  and  surveillance 
missions,  this  requirement  is  expected  to  grow  to  100  kilowatts,  (_1).  In 
addition,  a  space  nuclear  power  system  concept  with  approximately  100 
kilowatts  has  been  defined  for  such  uses  as  planetary  exploration, 
directed  energy  capabilities,  surveillance,  command,  control  and  com¬ 
munication  satellites,  etc.  Strategic  Defense  Initiative  (SDI)  programs 
have  also  identified  a  wide  range  of  power  requirements  with  standby 
power  needs  of  greater  than  10  kilowatts  and  peak  power  loads  greater 
than  10  megawatts.  In  this  context,  a  large  array  of  energy  sources  and 
power  conversion  concepts  (thermoelectric,  thermionic,  photo-voltaic, 
dynamic  isotope  power  systems,  etc.)  have  been  developed  and  proposed. 
The  present  thermal  control  methods  are  based  on  transport  of  heat  by 
solid  conductors  and  internal  radiation,  heat  pipes,  and  single-phase 
liquid  and  gaseous  loops.  The  high  operating  power  for  future  space 
applications  requires  more  efficient  thermal  transport  techniques.  Two- 
phase  loops  have  been  suggested  and  laboratory  tested  for  possible 
application  in  the  above  areas,  Refs.  (2)  and  (3).  In  comparison  to  a 
single-phase  loop,  the  two-phase  system  operates  at  considerably  smaller 
flow  rates  and  maintains  tighter  temperature  control  with  a  higher  heat 
transfer  coefficient.  However,  the  two-phase  fluid  flow  regimes, 
pressure  gradients,  and  heat  transfer  coefficients  must  be  evaluated  for 
application  in  the  weightless  environment  (zero-gravity)  of  an  orbiting 
satellite.  This  report  addresses  several  aspects  of  expected  flow  con- 
di ti ons . 


1.1  BACKGROUND 


Design  of  two-phase  heat  transport  systems  for  space  applications 
requires  a  knowledge  of  heat  and  mass  transfer  under  reduced  gravity 
conditions.  Many  of  the  issues  which  affect  system  performance  in  space 
have  been  identified  by  various  investigators.  These  include  two-phase 
flow  regime,  pressure  drop,  void-quality  relation,  and  two-phase  heat 
transfer.  In  some  instances,  the  technology  issues  identified  in  these 
studies  have  a  detrimental  effect  on  system  performance  and  also  pose 
severe  restrictions  on  system  efficiency. 

A  brief  look  at  the  simplified  forms  of  the  two-phase  flow  governing 
equations  will  reveal  the  need  and  importance  of  these  parameters. 


Void-Quality  Relationship: 


Generally  the  void-quality  relation  is  needed  to  reduce  the  number  of 
dependent  variables  in  the  two-phase  flow  conservation  equations.  A 
simplified  form  of  the  momentum  equation  in  the  absence  of  gravity  is 
given  below: 
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The  first  term  in  the  above  equation  is  the  friction  component  and  the 
second  term  is  the  inertia  or  acceleration  component.  The  flow  quality, 
x,  can  usually  be  determined  from  the  heat  input  to  the  system,  but  a 
void-quality  relation  is  needed  to  solve  for  the  inertia  and  therefore 
the  overall  pressure  loss. 


Two-Phase  Friction  Mu1tij>1  ier: 


Friction  loss  for  a  two-phase  system  is  generally  evaluated  by  using  the 
single-phase  relations  along  with  a  two-phase  friction  multiplier.  The 
relations  for  the  two-phase  friction  multiplier  at  reduced  gravities  are 
therefore  needed  for  design  of  the  spacecraft  cooling  systems. 

The  two-phase  friction  pressure  drop  can  be  expressed  in  terms  of  the 
single-phase  pressure  drop  as  follows: 
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<t>f0  is  the  two-phase  friction  multiplier  and  (dP/dx)  sp  is  the  singtle- 
phase  friction  pressure  drop  when  the  entire  flow  consists  of  liquid 
alone: 
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where  f  is  the  single-phase  friction  coefficient  and  D  is  the  channel 
diameter. 

Theoretically,  the  friction  multiplier  is  a  function  of  flow  quality, 
system  pressure,  and  mass  flux.  Therefore,  gravity  should  not  have  an 
effect  on  the  value  of  the  two-phase  multiplier.  However,  the  models 
developed  for  predicting  the  friction  multiplier  are  mostly  empirical 
and  it  is  not  clear  if  they  will  be  valid  in  the  absence  of  gravity  when 
the  phase  slip  and  distribution  of  voids  are  expected  to  be  different. 


Forced  Convective  Heat  Transfer  Coefficient: 


Prediction  of  the  thermal  performance  of  a  boiling/condensing  two-phase 
cycle  is  to  a  large  extent  dependent  on  the  accuracy  of  determining  the 
heat  transfer  coefficients.  The  most  desirable  mode  of  operation  for  a 
boiling  loop  is  the  forced  convective  nucleate  boiling  regime  where  the 
highest  heat  transfer  coefficient  is  obtained.  Among  the  existing  corre¬ 
lations,  the  most  successful  nucleate  boiling  heat  transfer  coefficients 
are  based  on  combining  a  forced  convective  single-phase  and  a  nucleate 
pool  boiling  heat  transfer  coefficient.  Empirical  factors  are  then 
introduced  to  account  for  the  influence  of  these  two  modes  of  heat 
transfer  under  different  system  conditions.  Applicability  of  these 
correlations  to  reduced  gravity  conditions  should  be  investigated. 

Two-Phase  How  Regime  Map: 

One  of  the  most  important  criteria  in  two-phase  flow  is  the  transition 
between  different  flow  regimes.  The  flow  regime  map  and  the  transition 
criteria  are  needed  in  order  to  operate  the  system  under  the  lowest  loss 
conditions  and  to  avoid  flow  regimes  which  result  in  instabilities  (slug 
flow).  In  addition,  the  two-phase  flow  regime  also  dictates  the  mode  of 
heat  transfer  which  will  determine  the  thermal  performance  of  the 
system. 

Current  knowledge  of  two-phase  flow  and  heat  transfer  is  mainly  derived 
from  terrestrial  experiments.  Unlike  pool  boiling  which  has  been  studied 
extensively  under  reduced  gravity  conditions,  very  little  work  has  been 
done  on  understanding  and  modeling  two-phase  flow.  One  of  the  reasons 
for  the  lack  of  information  at  microgravity  is  the  difficulty  in  per¬ 
forming  such  tests.  A  majority  of  low  gravity  fluid  mechanics  and  heat 
transfer  experiments  have  been  performed  in  drop  towers  or  under  simu¬ 
lated  reduced  gravity  conditions.  As  will  be  discussed  later,  studies 
have  been  recently  initiated  by  the  Air  Force,  NASA,  and  DOE  to  better 
understand  the  two-phase  flow  behavior  under  reduced  gravity  conditions. 


These  include  experiments  which  simulate  reduced  gravity  two-phase 
distribution,  drop  tower  tests,  and  the  airplane  trajectory  tests. 

Since  the  field  of  two-phase  flow  in  reduced  gravity  is  fairly  new,  most 
of  the  research  programs  have  aimed  at  developing  qualitative  infor¬ 
mation  mainly  for  confirming  the  existence  of  the  phenomena.  The  experi¬ 
mental  results  are  mainly  from  the  ground-based  tests  which  have  very 
short  duration.  Flow  development  and  existence  of  steady  state  con¬ 
ditions  for  such  a  short  duration  is  questionable.  Modeling  efforts  are 
limited  to  determination  of  the  flow  regime  transition  criteria.  The 
pressure  drop,  void  distribution,  and  heat  transfer  characteristics  have 
not  been  addressed  in  the  existing  studies. 

1.2  BRIEF  REVIEW  OF  PREVIOUS  EFFORTS 

The  subject  of  two-phase  flow  and  heat  transfer  has  been  studied  exten¬ 
sively  due  to  its  application  in  chemical  and  nuclear  industries.  It 
has  a  history  of  over  four  decades  and  there  are  numerous  text  books 
which  have  covered  different  aspects  of  two-phase  flow.  A  recent  book  on 
this  subject.  Ref.  (£),  contains  the  majority  of  useful  data  and 
modeling  efforts.  Due  to  the  complex  nature  of  two-phase  flow,  most  of 
the  problems  are  treated  empirically  and  even  the  mechanistic  approaches 
are  augumented  by  empirical  factors.  The  empirical  treatment  of  two- 
phase  flow  makes  extrapolation  to  other  conditions  e.g.,  reduced  gra¬ 
vity,  impractical. 

The  subjects  of  bubble  formation  and  pool  boiling  under  reduced  gravity 
conditions  have  been  studied  extensively.  Most  of  these  studies  have 
been  experimental  investigations  which  were  carried  out  in  drop  towers, 
(5)  to  (£).  Some  tests  have  been  performed  under  simulated  reduced  gra¬ 
vity  conditions  which  were  achieved  by  compensating  for  the  earth's  gra¬ 
vity  by  using  strong  magnetic  force,  Refs.  (£0)  and  (£1).  Basically  the 
pool  boiling  experiments  have  shown  negligible  change  in  the  nucleate 


boiling  heat  transfer  coefficient  and  considerable  decrease  in  critical 
heat  flux  with  decreasing  gravities.  Bubble  formation  and  departure  has 
been  studied  by  resolution  of  gravity  into  components  in  a  tilted  (or 
horizontal)  container.  Ref.  (12)  and  (_13).  The  theory  of  vanishing 
nucleate  boiling  in  the  absence  of  gravity.  Ref.  (_14),  has  yet  to  be 
proven  in  future  experiments. 

Unlike  the  areas  of  pool  boiling,  critical  heat  flux,  and  bubble  growth 
mechanism,  very  little  work  has  been  done  on  understanding  and  modeling 
two-phase  flow  at  reduced  gravities.  A  review  article  by  Siegel,  Ref. 
( 15) ,  has  covered  the  data  available  from  the  reduced  gravity  experi¬ 
ments  prior  to  1967.  A  summary  of  pre-1974  literature  can  also  be  found 
in  Ref.  (16).  The  most  recent  review  of  the  experimental  efforts  in 
reduced  gravity  two-phase  flow  can  be  found  in  Ref .(17). 

Results  of  an  experimental  study  on  pressure  and  temperature  changes  in 
forced  convective  boiling  at  zero  gravity  have  been  reported  in  Ref. 
( 18).  It  has  been  shown  that  the  system  pressure  and  temperature 
increase  and  boiling  oscillations  damp  out  in  zero  gravity.  Low  heat 
flux  convective  boiling  at  zero  gravity  using  a  drop  tower  facility  was 
studied  in  Ref.  (T9)  and,  mainly,  the  bubble  diameter  was  determined. 
Other  two-phase  flow  studies  have  reviewed  the  concept  of  two-phase  flow 
application  for  spacecraft,  (20)  and  (2J_).  Since  these  efforts  have 
basically  studied  the  qualitative  nature  of  two-phase  flow  at  microgra¬ 
vity,  it  is  anticipated  that  the  most  useful  information  will  be 
obtained  from  the  programs  which  are  presently  in  progress.  These 
efforts,  which  have  been  initiated  by  the  Air  Force,  NASA  and  DOE,  are 
expected  to  be  completed  in  late  1988.  A  brief  summary  of  the  objectives 
of  the  major  efforts  in  progress  is  given  below. 


AETA/Dartmouth/University  of  New  Hampshire  study.  Ref.  (22): 

This  is  an  experimental  and  analytical  study  of  two-phase  flow  regime 
transitions  under  reduced  gravity  conditions.  The  outcome  is  expected  to 
consist  of  experimental  data  for  simulated  reduced  gravity  conditions 
(equal  density  liquids  at  earth  gravity)  and  a  model  for  two-phase  flow 
regime  transition. 

NASA  *  DOE  Supported  University  Activities: 

NASA  and  DOE  have  contracts  with  three  universities  to  perform  a  series 
of  tests  with  facilities  ranging  from  drop  towers  to  Spacelab.  These 
efforts  are: 

a)  University  of  Houston  Study,  Ref.  ( 23) : 

This  is  an  experimental  and  analytical  study  of  two-phase  flow  regime 
transition,  interfacial  instabilities,  and  slug  flow  modeling.  The  tests 
are  performed  with  air-water  in  the  NASA  drop  tower  facility  and  aboard 
the  NASA  Learjet. 

b)  University  of  Michigan  Study: 

This  university  has  a  program  to  study  pool  boiling  (with  the  university 
drop  tower)  and  convective  boiling  (aboard  NASA  airplane).  This  study 
will  mostly  concentrate  on  low  quality  regions  including  incipience  of 
boiling.  Tests  are  planned  for  flat  surfaces  and  eventually  inside 
tubes . 

c)  Texas  A4M  University: 

This  is  an  experimental  program  to  study  forced  condensation  and  heat 
transfer  rates.  The  tests  are  carried  out  in  the  NASA  KC-135  aircraft. 

Renewed  interest  in  the  space-based  nuclear  reactor  power  systems  has 
initiated  new  areas  in  the  field  of  reactor  thermal -hydraul ics  and 
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therefore  two-phase  flow  and  heat  transfer.  National  laboratories  and 
private  industry  have  started  to  evaluate  the  two-phase  constitutive 
models  which  are  needed  in  the  computer  codes  used  for  steady  and  tran¬ 
sient  analysis  of  the  reactors,  ( 24 )  (25)  (_26).  The  Battelle  study,  Ref. 
(24) ,  includes  an  experimental  study  of  forced  convective  boiling  which 
was  carried  out  in  the  NASA  drop  tube  and  aboard  the  NASA  KC-135.  A 
total  of  39  aircraft  tests  were  completed  in  1986  and  the  data  report 
will  be  available  in  late  1987.  In  addition,  the  possibility  of  using 
the  existing  computer  codes  for  analysis  of  the  reactor  performance 
under  reduced  gravity  conditions  was  investigated. 

A  series  of  tests  with  convective  boiling/condensing  water,  and  under 
adiabatic  conditions  with  nitrogen/water  at  atmospheric  conditions  were 
supported  jointly  by  the  Air  Force  Wright  Aeronautical  Laboratories  and 
the  Johnson  Space  Center  aboard  the  NASA  KC-135,  Ref.  (27).  In  addition 
to  flow  visualization  which  provided  information  on  flow  regime  tran¬ 
sitions,  these  tests  provided  pressure  drop  data  for  low  flow  rates. 

NASA  JSC  has  the  responsibility  for  thermal  management  of  space  station 
equipment  external  to  the  pressurized  modules.  Under  this  program, 
several  prototype  development  activities  for  the  Thermal  Bus  System 
(TBS)  were  initiated.  The  Boeing  TBS  concept  uses  a  rotary  fluid  manage¬ 
ment  device  for  separation  of  liquid  from  the  two-phase  mixture  exiting 
from  the  evaporators.  The  Grumman  concept  controls  the  flow  into  evap¬ 
orators  to  obtain  a  high  quality  vapor  at  the  exit.  In  order  to  gain  a 
better  understanding  of  the  two-phase  flow  phenomena  to  ensure  adequate 
and  efficient  hardware  design,  several  ground-based  studies  including 
the  university  activities  mentioned  above  and  a  KC-135  testing  of  a  two- 
phase  loop  developed  by  Sundstrand  were  initiated  by  NASA.  This  loop 
will  be  used  for  identification  of  the  two-phase  flow  regimes  and 
measurement  of  pressure  drops. 

Design  of  a  facility  for  performing  a  series  of  two-phase  flow  experi¬ 
ments  aboard  the  shuttle  vehicle  has  been  under  way  for  over  eight 
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years.  Refs.  (28)  (29).  A  separate  study  for  preliminary  design  of  a 
similar  two-phase  experiment  has  been  initiated  by  NASA  JSC,  (30) . 

1.3  OBJECTIVES  AND  SCOPE  OF  EFFORT 

This  project  was  funded  by  the  Air  Force  Wright  Aeronautical 
Laboratories  as  a  Phase  II  SBIR  program  to  study  two-phase  flow  behavior 
under  zero  gravity  conditions.  The  overall  objectives  of  this  study  were 
to  generate  a  data  base  for  two-phase  pressure  drop  and  void-quality 
relation  under  simulated  zero  gravity  conditions  and  to  develop  analyti¬ 
cal  models  to  predict  these  parameters  for  bubbly  and  annular  flow. 

The  simulation  of  zero  gravity  two-phase  flow  was  achieved  by  using  two 
immiscible  liquids  with  equal  densities  to  eliminate  the  buoyancy  com¬ 
ponent.  Although  this  approach  does  not  eliminate  the  gravity  effects, 
it  does  provide  a  representation  for  void  distribution  in  the  absence  of 
gravity.  Water  was  used  as  one  of  the  working  fluids  and  simulated  the 
gaseous  phase,  while  the  second  fluid  was  selected  to  simulate  the 
liquid  phase.  The  property  ratios  of  the  two  liquids  were  in  the  range 
of  the  two-phase  liquid/gas  flows.  The  tests  were  performed  under  steady 
state  fully  developed  turbulent  flow  conditions,  since  it  is  believed 
that  for  most  practical  applications  the  two-phase  heat  transfer  loops 
operate  within  the  turbulent  flow  range.  The  major  parameters  are 
pressure  drop  and  cross  sectional  average  void  fraction,  and  the  test 
variables  consist  of  total  flow  rate,  quality,  and  fluid  properties. 

The  modeling  effort  is  limited  to  developing  relations  for  the  two-phase 
friction  multiplier  and  void-quality  relation  under  bubbly  and  annular 
flow  conditions.  The  bubbly  flow  model  is  based  on  the  assumption  of 
local  homogeneous  conditions  between  the  phases  but  allows  for  void 
distribution  in  the  radial  direction.  Separated  flow  conservation 
equations  are  used,  and  single-phase  turbulent  flow  eddy  diffusivity 
relations  are  employed.  In  addition,  a  single  particle  lift  force  model 


is  utilized  to  account  for  radial  forces  on  the  bubbles  within  the  shear 
layer.  A  tentative  model  is  proposed  for  bubbly  flow  with  nucleation  at 
the  wall.  This  approach  uses  the  model  developed  for  adiabatic  flow  but 
the  boundary  condition  is  modified  to  account  for  the  presence  of  voids 
at  the  wall.  The  annular  flow  model  is  based  on  the  triangular  approach 
of  Hewitt,  Ref.  (31.)*  The  interfacial  shear  and  entrainment  relations 
are  empirically  determined  from  the  results  of  the  experimental  effort 
with  equal  density  liquids. 

It  should  be  noted  that  in  the  present  analytical  effort  no  attempt  was 
made  to  model  the  heat  transfer  characteristics  of  two-phase  flow  in 
zero  gravity.  It  is  believed  that  the  mechanism  of  heat  transfer  will 
not  be  affected  by  the  absence  of  gravity.  However,  the  heat  transfer 
rates  are  influenced  by  void  distribution  within  the  two-phase  flow.  The 
earth  gravity  models  are  generally  based  on  flow  quality  and  therefore 
cannot  be  applied  to  reduced  gravity  conditions.  If  one  can  incorporate 
the  differences  in  void-quality  relation  into  the  earth  gravity  correla¬ 
tions,  they  can  be  used  under  reduced  gravity  conditions. 


Section  2 

EXPERIMENTAL  INVESTIGATIONS 

The  experimental  program  was  aimed  at  investigating  the  two-phase  flow 
behavior  under  zero  gravity  conditions.  The  tests  were  planned  to  simu¬ 
late  the  phase  distribution  and  two-phase  pressure  drop  characteristics 
of  gas-liquid  flow  in  the  absence  of  gravity.  The  major  objective  of 
this  program  was  to  generate  a  data  base  for  the  two-phase  friction 
multiplier  and  void-quality  relation  under  simulated  bubbly  and  annular 
flow  patterns.  This  was  accomplished  with  a  flow  of  two  immiscible 
liquids  of  nearly  equal  densities.  Although,  equal  density  liquid-liquid 
flow  will  not  eliminate  the  gravity  effect,  it  should,  however,  be  a 
close  representation  of  the  phase  distribution.  Proper  scaling  would 
also  require  a  simulation  of  the  viscosity  and  interfacial  tension 
ratios. 

The  major  test  parameters  were  the  void  fraction  (or  hold-up  ratio)  and 
the  pressure  drop  along  the  test  section.  The  test  variables  included  the 
total  mass  flow  rate,  relative  phase  content  (quality),  and  the  fluid 
properties. 

2.1  SELECTION  OF  LIQUIDS 

Distilled  water  was  chosen  to  be  one  of  the  test  fluids.  Because  of  the 
low  viscosity  of  water  compared  with  other  liquids,  it  was  convenient  to 
have  water  simulate  the  gas  in  a  zero-gravity  gas-liquid  flow  simula¬ 
tion.  For  the  second  fluid  to  properly  simulate  a  liquid  relative  to 
water,  a  suitable  viscosity  for  the  fluid  must  be  considered.  The 
distilled  water  used  was  of  density  0.997  g/ml  and  kinematic  viscosity 
0.894  cSt  at  atmospheric  pressure  and  room  temperature  (25°C).  The  ratio 
of  kinematic  viscosity  for  water  and  air  at  the  above-mentioned  con¬ 
ditions  is  about  15.  Based  on  this  information,  Dimethyl  silicone  fluid 
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L-45,  Union  Carbide  Corporation  and  Dowtherm  LF  heat  transfer  fluid 
(DT-LF)  were  chosen  to  simulate  the  liquid  respectively  in  different 
tests.  The  properties  of  the  test  fluids  are  listed  below: 

Fluid  Properties  at  Room  Temperature 


VISCOSITY 

SPECIFIC 

SURFACE  TENSION 

RESISTIVITY 

cSt 

GRAVITY 

dyne/cm 

ohm-cm 

WATER 

0.894 

1.0 

71 .95(water-ai r ) 

2xl04 

L-45 

10.0 

0.938 

32.31 (water-L-45) 

lxlO14 

DT-LF 

4.323 

1.038 

48.57(water-DT-LF) 

6xlOH 

Since  turbulent  forces  dominated  the  buoyancy  effects  in  this  test, 
specific  gravity  values  of  0.938  and  1.038  were  considered  to  be  suf¬ 
ficiently  close  to  one  for  simulating  the  two-phase  flow  character 
under  zero  gravity  condition. 

2.2  TEST  LOOP  AND  INSTRUMENTATION 

The  experimental  apparatus  used  to  simulate  horizontal  gas-liquid  pipe 
flow  under  zero  gravity  conditions  is  shown  in  Figure  2.1 

Distilled  water  (white  arrow)  and  simulated  liquid  (L-45  or  DT-LF) 
(black  arrow)  were  pumped  from  separate  storage  tanks  into  the  mixing 
section.  Different  mixing  sections  were  designed  to  achieve  the  desired 
flow  regimes.  The  flow  of  the  two  fluids  (black-and-white  arrow)  simu¬ 
lating  two-phase  flow  was  produced  in  the  mixing  section  and  flowed  into 
a  twelve-foot-long  horizontal  test  section.  Pressure  drop  and  void 
fraction  measurement  units  were  installed  in  this  section.  After  the 
test  section,  the  two-phase  mixture  flowed  into  a  receiving  tank  where 
it  was  separated. 
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2.2.1  Tanks 

The  storage  tanks  for  the  distilled  water  and  the  simulated  liquid  are 
cylindrical  graduated  HDPE  tanks  of  volume  150  gallons.  All  tanks  were 
from  Cole-Parmer  Instrument  Company. 

2.2.2  Pumps 

The  water  pump  is  a  PRICE  1.5  hp  centrifugal  plastic  impeller,  bronze 
body  pump  (PRICE  model  HP  75-100)  with  a  maximum  capacity  of  18  gpn  at 
120  feet  head. 

2.2.3  Rotameters 

The  two  rotameters  are  Fischer  and  Porter  glass  tube  indicating  rotame¬ 
ters.  The  rotameter  used  for  simulated  liquid  had  a  viscosity  limit  of 
18.5  cp  and  a  maximum  flow  rate  of  31.6  gpn.  The  rotameter  used  for 
water  had  a  maximum  capacity  of  17  gpn  at  70°F. 

2.2.4  Mixing  Sections 

In  order  to  simulate  both  bubbly  flow  and  annular  flow,  two  different 
sections  were  made  to  meet  the  purpose. 

The  bubbly  flow  mixing  section  is  shown  in  Figure  2.2.  Distilled  water 
was  introduced  through  120  1/64  inch  diameter  holes  on  a  perforated  tube 
which  was  centrally  located  in  the  interior  of  simulated  liquid  flow. 
The  distilled  water  emerging  from  the  injection  holeswas  carried  away 
from  the  tube  wall  by  the  simulated  liquid  flow,  and  water  bubbles  were 
formed . 

The  annular  flow  mixing  section  is  shown  in  Figure  2.3.  Simulated 
liquid  was  pumped  into  an  annular  area  of  annulus  thickness  0.075  inch 
to  form  a  simulated  liquid  film  on  the  inner  wall  of  the  test  section. 
The  distilled  water  (simulating  the  gas)  was  introduced  into  the  core, and 


the  result  was  an  annular  type  flow,  i.e.  simulated  liquid  on  the  wall 
and  gas  (simulated  by  water)  in  the  core.  This  test  section  was 
designed  to  make  the  film  thickness  interchangeable.  The  mixing  section 
for  L-45  water  flow  was  made  of  plexiglass,  the  one  for  DT-LF  water  flow 
was  made  of  pvc  pipes. 

2.2.5  Test  Section 

The  test  section  was  a  twel  ve-foot-1  ong ,  one  inch  outer  diameter,  3/4 
inch  inner  diameter  round  tube  of  which  the  axis  was  in  the  horizontal 
direction.  For  L-45  water  flow,  clear  plexiglass  tube  was  used,  for 
DT-LF  water  flow,  pc  pipe  was  used.  Pressure  drop  and  void  fraction 
measurements  were  carried  out  on  this  section.  The  test  section  was 
checked  to  be  hydrodynamical ly  smooth  by  a  single-phase  friction  factor 
measurement. 

2.2.6  Pressure  Drop  Measurement  Unit 


The  pressure  drop  measurement  unit  consists  of: 

(a)  five  pressure  taps 

(b)  control  panel 

(c)  two  pressure  transducers 

(d)  digital  transducer  indicator 

Five  2x2x1  inch  plexiglass  blocks  were  equally  spaced  along  the  test 
section  to  provide  for  the  pressure  taps. 

These  were  spaced  at  a  distance  two  feet  and  nine  inches  between  two 
neighboring  taps.  1/8  inch  diameter  holes  were  drilled  through  the  wall 
of  the  test  section  and  the  plexiglass  blocks  to  record  the  pressure. 
The  five  pressure  taps  were  connected  to  five  globe  valves  on  the 
control  panel  to  control  the  measurements  of  the  pressure  drops  which 
were  measured  between  pressure  tap  number  one  and  the  other  four  taps. 
The  configuration  of  the  pressure  taps  is  shown  in  Figure  2.4.  The  cir¬ 
cuits  of  the  control  panel  are  shown  in  Figure  2.5. 


POT 


Figure  2.5  -  Pressure  Sensing  Unit 


Two  Validyne  DP15  variable  reluctance  differential  pressure  transducers 
were  used  to  measure  the  pressure  difference  between  two  pressure  taps. 
The  two  transducers  were  equipped  with  appropriate  diaphrams  for  maximum 
pressure  ranges  from  zero  to  0.8  psid  and  from  zero  to  8.0  psid,  respec¬ 
tively.  The  readings  of  the  measured  results  were  obtained  from  a 
Validyne  model  CD  223  dual  channel  digital  transducer.  A  switch  allowed 
toggling  between  the  readings  from  the  two  transducers. 

2.2.7  Void  Fraction  Measurement  Unit 

The  void  fraction  (hold-up  ratio)  measurement  unit  consists  of: 

(a)  two  quick-closing  valves  and  actuators 

(b)  two  actuation  gas  tanks 

(c)  one  nitrogen  supply  tank 

(d)  impedance  probe 

Two  Worcester  Controls  Series  59,  one  inch  full  port  ball  valves  with 
brass  body  and  teflon  seal  are  spaced  a  distance  of  twenty  inches  bet¬ 
ween  pressure  taps  number  three  and  number  four  for  trapping  the  two- 
phase  mixture  between  them.  The  volume  of  the  fluid  trapped  by  the 
valve  was  150  ml.  The  two  valves  were  actuated  simultaneously  by  two 
Kinetrol  model  05  vane  type  double  acting  air  actuators  with  closing 
time  less  than  one  second.  The  two  actuators  were  supplied  with  nitro¬ 
gen  by  the  two  actuation  gas  tanks  which  were  charged  with  a  high 
pressure  nitrogen  supply  tank.  The  trapped  mixture  could  be  drained  out 
through  a  half-inch  drain  by  high  pressure  nitrogen  for  the  void  frac¬ 
tion  measurement. 

An  Auburn  International  Capacitance  gauge  model  1090  was  used  to  provide 
an  alternative  determination  of  the  void  fraction  by  measuring  the 
variation  of  the  capacitance  of  the  flow.  The  detailed  working  prin¬ 
ciples  of  this  capacitance  probe  are  given  in  the  Appendix. 
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2.3  TEST  MATRIX 

The  two-phase  flow  condition  in  the  test  was  specified  in  terms  of  the 
total  Reynolds  number  (Ret,  based  on  simulated  liquid  properties)  and 
the  mean  phase  content  or  mass  quality  (x,  water  mass  flow  rate  divided 
by  the  total  mass  flow  rate).  The  selected  total  Reynolds  numbers  and 
mass  qualities  for  L-45/water  flow  and  DT-LF/water  flow  with  bubbly  and 
annular  flow  mixing  sections  are  listed  below: 

Bubbly  Flow  Mixing  Section 


Ret 

Min(x) 

L-45/Water  Flow 

Max(x) 

#  of 

data 

5,500 

0.13 

0.80 

13 

7,900 

0.08 

0.68 

20 

8,800 

0.11 

0.51 

7 

10,200 

0.07 

0.45 

5 

11,600 

0.11 

0.37 

5 

Ret 

Mi  n  ( x ) 

DT-LF/Water  Flow 

Max(x) 

#  of 

data 

9,000 

0.17 

0.45 

5 

15,900 

0.09 

0.41 

5 

18,000 

0.08 

0.36 

5 

8 
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Annular  Flow  Mixing  Section 


L-45/Water  Flow 

Annulus  Thickness  =  0.191  cm 
Ret  Min(x)  Max(x)  #  of  data  points 

6.900  0.39  0.84  6 

7.900  0.38  0.74  5 

8.900  0.40  0.66  5 

Annulus  Thickness  =  0.095  cm 
Ret  Min(x)  Max(x) 

5,500  0.54  0.79 

2.4  TEST  PROCEDURE 

Before  the  test,  the  two  actuation  gas  tanks  in  the  void  fraction 
measurement  system  were  charged  with  nitrogen  to  sixty  psig  for  later 
use.  The  impedance  probe  was  calibrated  by  running  single-phase 
distilled  water  and  simulated  liquid  respectively  through  it  to  set  the 
zero  value  (when  simulated  liquid  was  in  the  probe)  and  the  full  scale 
value  (when  measuring  only  the  distilled  water). 

For  every  given  total  mass  flow  rate,  the  simulation  started  with  running 
simulated  liquid  through  the  whole  test  section  to  build  up  the  "liquid 
phase"  matrix.  Distilled  water  was  then  introduced  into  the  simulated 
liquid  flow  through  the  mixing  section  to  produce  the  desired  flow  pat¬ 
tern.  The  proper  mass  quality  was  obtained  by  adjusting  the  control 
valves  for  the  two  fluids.  The  pressure  drops  between  the  first  pressure 
tap  and  the  other  four  taps  were  measured  with  the  pressure  transducer 
system  discussed  in  Section  2.2.6.  The  measurement  of  the  local  void 
fraction  between  the  third  and  fourth  pressure  taps  could  be  carried  out 


#  of  data  points 
5 


2-12 


in  two  ways.  One  was  by  the  impedance  probe  described  in  2.2.7 ,  and  the 
other  method  was  to  measure  the  void  fraction  directly  by  removing  some 
of  the  mixture  from  the  test  section.  To  accomplish  this,  the  two  pumps 
were  shut  off  to  stop  the  flow  of  the  two  fluids  and  right  after  this, 
the  solenoids  on  the  quick-closing  valves  were  activated  to  close  the 
valves.  This  contained  the  water-simulated  liquid  mixture  between  the 
two  valves.  This  also  terminated  the  test  and  the  isolated  mixture  was 
drained  for  the  void  fraction  measurement. 

To  completely  drain  the  fluid  between  the  quick-closing  valves  into  a 
graduated  cylinder,  the  drain  was  open  and  nitrogen  of  forty  psig  was 
introduced  into  the  isolated  section  to  drain  the  fluid  into  a  graduated 
cyl i nder. 

After  the  mixture  sample  had  been  drained,  the  system  was  ready  for 
another  test.  The  temperature  of  the  flow  was  determined  by  measuring 
the  temperature  of  the  mixture  in  the  receiving  tank  after  the  pumps  had 
been  shut  off. 

In  order  to  separate  the  simulated  liquid  from  the  distilled  water  in 
the  graduated  cylinders  and  the  receiving  tank,  three  to  four  days  were 
needed  for  the  two  liquids  to  be  separated  by  gravity.  After  the 
separation  had  been  accomplished,  the  distilled  water  was  drained  and 
discarded.  The  simulated  liquid  was  pumped  back  through  the  bypass  line 
into  its  storage  tank  for  future  use. 

2.5  DATA  REDUCTION 

The  test  results  are  presented  in  the  form  of  two-phase  frictional 
multiplier  (<tf2),  based  on  the  liquid  phase,  namely  L-45  or  DT-LF,  and 
void  fraction  (a)  as  functions  of  Martinelli  parameter  (X)  respectively, 
to  describe  the  relations  between  pressure  losses  and  flow  conditions. 

In  order  to  relate  these  two-phase  variables  to  the  measured  quantities, 
the  following  calculations  were  performed: 
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Two-phase  Variables 


mass  quality 

two-phase  frictional  multiplier 
Martinelli  parameter 


Measured  Quantities 


ql  :  simulated  liquid  flow  rate  (gpm) 
qw  :  distilled  water  flow  rate  (gpm) 
p2,p3,p4,p5  :  pressure  drops  between  pressure  tap 
number  1  and  pressure  taps  number  2, 3, 4, 5  (psi) 


Fluid  Properties- 


pw  =  0.997  g/cn3  Density  of  Water 

p$  =  0.935  g/cm3  Density  of  1-45 

Pg  =  1.035  g/cm3  Density  of  DT-LF 

Uy  =  0.891  cPs  Viscosity  of  Water 

us  =  93.50  cPs  Viscosity  of  L-45 

UD  =  44*75  cPs  Viscosity  of  DT-LF 
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where 

dP 

( — )  :  two-phase  frictional  pressure  gradient 
d  1  f 


A  :  frictional  pressure  gradient  for  liquid 
dz  f 

phase 

(d£)  _  pressure  drop  between  taps  _  _AP 

dz  p  distance  between  taps  L 

In  order  to  minimize  the  entrance  effect  between  pressure  taps  number  1 
and  2  and  the  disturbance  by  the  capacitance  probe  between  pressure  taps 
number  3  and  4,  P  was  chosen  as  P  =  (P3  -  P2)  +  (P5  -  P4),  and  L  was 
taken  to  be  twice  the  distance  between  the  neighboring  pressure  taps, 
which  was  66  inches. 

Therefore 


-  (IP)  ,  P3.,t_P5JL.P2.-.P4  (psi/irK) 
dz  p  66 

The  single-phase  frictional  pressure  gradient  was  written  in  terms  of 
friction  factor  as: 


-  <£>  --  4  f 


dz'f"  0  ’ 


2 


where 

D  =  the  diameter  of  the  pipe  (1.905  cm) 

^ i  =  the  single-phase  friction  f^'tor  of  the  simulated  liquid 

u  i  =  the  mean  velocity  of  the  simulated  liquid  when  it  flows  alone 
in  the  pipe 
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By  the  relation 


q,  _  4q, 
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A  ttD 


k 

3 
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we  have 


-  (f ) 

dz  f 


32  f  „  2 
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2.5.3  Relation  Between  and  Measured  Quantities 


Definition 


dP 

^dzV 

W 

02  n 


Where 

=  frictional  pressure  gradient  for  gas  phase  alone 

dZ  g 


ft 

dz  f 


tptuz2 


With 


/dP  %  _  32  f  n  2 

dz  f  ,,0 


Similarly 


Ww 


2 


where  f  is  the  single-phase  friction  factor  for  water 

W 


dP 

(dz)f 


=  (_*)  (JL)  (-L) 


When  L-45  is  used: 


/fs  % 
Xtt  =  0.968  / —  — 

u  V  f  -  q 

>  w  Mw 


When  DT-LF  is  used: 

'  w  Mw 

2.5.4  Calculation  of  Single-Phase  Friction  Factors,  fj i,  fw 


The  single-phase  friction  factor  for  pipe  flow  is  usually  written  as  a 
function  of  Reynolds  number  Re,  defined  as  Re  =  puD/u,  for  both  laminar 
and  fully  developed  turbulent  flow.  Because  the  range  of  the  Reynolds 
numbers  for  simulated  liquid  and  distilled  water  covers  the  transition 
region  (for  smooth  pipe  2,000  <  Re  <  8,000),  and  there  is  no  single 
relation  between  friction  factor  and  Reynolds  number  suitable  to  be  used 
for  the  calculation  of  f,  the  experimental  results  (Figure  2.6)  were 
used  for  this  Reynolds  number  range.  The  experimental  data  were  divided 
into  five  sections:  Re  <  2,000,  2,000  <  Re  <  4,000,  4,000  <  Re  <  6,800, 
6,800  <  Re  <  8,000,  and  8,000  <  Re.  The  results  in  the  middle  three 
sections,  which  covered  the  laminar-turbulent  transition  region,  were 
approximated  by  cubic  polynomials  obtained  by  the  difference  table 
method.  The  results  are  given  as: 


/.v.v.v.v.v.v;. 


Figure  2.6  -  Friction  Factor  for  Pipe  Flow  in  Laminar  - 
Turbulent  Transition  Region,  from  Compact 
Heat  Exchangers,  W.  L.  Kays  and  A.  L.  London 
McGraw-HilK 
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Re  <  2,000:  } 

f  -ii  : 

Re  ; 

2,000  <  Re  <  4,000:  j 

\ 

X  -  Re  ' 

”  1,000 

X2  =  X-2  j 

X3  =  X-2. 5  ■ 

- 

f  =  10"5  x  (902-304{ X2)  +  148(X2)(X3)  -  34.35( X2) ( X3) ( X-3) ) 

4,000  <  Re  <  6,800: 

X  -  Re 

"  1,000 

X4  =  X-4 

X5  =  X-5 

X6  =  X-6 

v  =  10"5  (635  +  90 ( X4)  -  22.5(X4)(X5)  +  2.083( X4) ( X5) ( X6) ) 

6,800  <  Re  <  8,000: 

X  =  Re 

1 ,000 

X6  =  X-6. 8 

X7  =  X-7 

X8  =  X-7. 5 


f  =  10-5 


(782  +  40(X6)  -  100(X6)(X7)  +  83.33( X6) X7) X8) ) 


isriv 


rLi*. 


8,000  <  Re: 


f  =  0.05858  x  Re"0,2258 

To  relate  the  single-phase  Reynolds  numbers  of  the  test  fluids  to  the 
measured  quantities,  start  with  the  definition 


Since 


q  =  uA  =  -  D2u 
4 


Re  =  lea 

■nuD 


1  gpm  =  63.09  cm  /sec 
D  =  0.75  in.  =  1.905  cm 


We  have  for  water  flow 
Rew  =  4718.35  q^ 


For  L-45  flow 

Res  =  421.67  qc 


For  DT-LF  flow 

ReD  =  975.44 
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2.6  TEST  RESULTS 


2.6.1  Two-Phase  Multiplier 


The  two-phase  multipliers  for  different  flow  conditions  are  shown  in  the 
form  of  two-phase  multiplier  versus  Martinelli's  parameter  in  Figures 
2.7  through  2.23  and  are  summarized  in  Table  2.1 

A  correlation  for  the  two-phase  multiplier  as  a  function  of  Martinelli's 
parameter  based  on  all  the  experimental  data  (Figure  2.23)  was  proposed 
as  follows: 


l 


(2.1) 


This  relation  is  plotted  as  solid  curves  on  the  data  figures.  For  com¬ 
parison,  Martinelli's  correlation  for  air-water  flow  (33): 


1  +  20 _  +  J_ 

I  Y  v  2 
tt  Xft 


(2.2) 


is  also  presented  on  the  figures  as  dashed-dotted  curves. 

2.6.2  Void  Fraction 

The  measured  void  fractions  are  plotted  versus  Martinelli's  parameter 
for  different  flow  conditions  in  Figures  2.24  to  2.36  and  are  summarized 
in  Table  2.2.  The  recommended  correlation  for  void  fraction  as  a  func¬ 
tion  of  inverse  Martinelli's  parameter  based  on  all  the  experimental 
data  (Figure  2.36)  is  obtained  as: 


i 


« -  (1  ♦  (/-l'1'45)'0-8 

xtt 


(2.3) 


and  is  plotted  on  the  figures  in  solid  curves.  Martinelli's  correlation 
for  void  fraction  based  on  air-water  flow  (33)  is  written  as: 


a 


(i  *  a.)-0-8)"0-378 

*tt 


(2.4) 


which  is  shown  as  dashed-dotted  curves  on  the  figures. 
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Table  2.1 


FLOW  CONDITIONS  FOR  TWO-PHASE  MULTIPLIER  FIGURES 


igure 

Ret 

Test  Fluids 

Mixing  Section 

2.7 

All 

L-45-water 

Bubbly  and  Annular 

2.8 

All 

L-45-water 

Bubbly 

2.9 

All 

L-45-water 

Annul ar 

2.10 

5,500 

L-45-water 

Bubbly 

2.11 

7,900 

L-45-water 

Bubbly 

2.12 

8,800 

L-45-water 

Bubbly 

2.13 

10,200 

L-45-water 

Bubbly 

2.14 

11,600 

L-45-water 

Bubbly 

2.15 

6,900 

L-45-water 

Annul ar 

2.16 

7,900 

L-45-water 

Annul ar 

2.17 

8,900 

L-45-water 

Annul ar 

2.18 

5,500 

L-45-water 

Annular 

2.19 

All 

DT -LF -water 

Bubbly 

2.20 

9,000 

DT-LF -water 

Bubbly 

2.21 

15,900 

DT-LF -water 

Bubbly 

2.22 

18,000 

DT-LF-water 

Bubbly 

2.23 

All  the  data  points 
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Test  Fluids 


Mixing  Section 


L-45-water 
L-45-water 
L-45-water 
L-45-water 
L-45-water 
L-45-water 
L-45-water 
L-45-water 
DT-LF -water 
DT-LF -water 
DT-LF -water 
DT -LF -water 
All  the  data  points 


Bubbly  and  Annular 
Bubbly 
Annular 
Bubbly 
Bubbly 
Bubbly 
Bubbly 
Annular 
Bubbly 
Bubbly 
Bubbly 
Bubbly 
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Figure  2.22  Two-Phase  Friction  Multiplier  Versus  Inverse  Martinelli  Parameter 
for  DT-LF/Water  Flow  with  Bubbly  Mixing  Section  for  Total  Reynolds 
Number  of  Re  =  18000. 

The  Solid  Curve  corresponds  to  Equation  2.1  and  the  Dashed-Dotted 
Curve  to  Equation  2.2. 
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Section  3 

BUBBLY  FLOW  V0I0  DISTRIBUTION  AND  FRICTION  MULTIPLIER 


In  the  absence  of  gravity,  the  distribution  of  bubbles  in  a  bubbly  mix¬ 
ture  within  a  pipe  is  controlled  by  the  drag  force,  the  lift  force,  and 
the  turbulent  fluctuations  in  the  liquid  phase.  Since  the  ratio  of  the 
surface  area  to  volume  of  the  bubbles  is  large,  the  bubbles  are  expected 
to  flow  at  nearly  the  same  velocity  as  the  liquid  in  the  axial  direc¬ 
tion.  However,  since  the  lift  force  and  turbulent  fluctuations  cause 
the  bubbles  to  migrate  towards  the  region  of  higher  velocity  (axis  of 
the  pipe),  a  void  distribution  will  be  established  in  the  radial  direc¬ 
tion.  Even  though  the  flow  can  be  assumed  to  be  in  a  local  homogeneous 
condition  (due  to  nearly  equal  phase  velocities),  homogeneous  flow 
assumptions  cannot  be  made  and  the  slip  ratio  would  be  different  from 
unity  due  to  the  non-homogeneous  void  distribution  in  the  radial  direc¬ 
ti  on . 

Separated  flow  conservation  equations  are  used  to  obtain  a  mathematical 
model  for  the  flow.  The  governing  equations  which  consist  of  the  momen¬ 
tum  conservation  for  the  liquid  and  vapor  phases,  and  mass  conservation 
for  the  liquid  phase  are  solved  to  develop  a  method  for  calculating  the 
void  distribution  and  two-phase  friction  multiplier  for  the  bubbly  flow 
under  zero  gravity  conditions.  The  approach  is  similar  to  the  method 
used  by  Wang,  Ref.  (34).  However,  the  lift  force  relations  used  are 
different,  and  in  the  present  study  several  assumptions  unique  to  the 
zero  gravity  conditions  can  be  made  which  enabled  us  to  develop  a 
general  solution  for  the  desired  parameters. 
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3.1  GOVERNING  EQUATIONS  AND  THE  METHOD  OF  ANALYSIS 


For  a  steady  adiabatic  two-phase  flow  with  negligible  interfacial  mass 
and  momentum  transfer,  the  general  form  of  the  conservation  equations  in 
a  two-dimensional  cartesian  coordinate  is: 


,(pkVk>  ,  3<°kakV 


(3.1) 


Momentum 


xdir.  2-  (»kVkZ>  <°k“kVk> 


y  dir.  f-  (rkVkuk>  ♦  §-  (»kvk2) 


3  P  9  /  , 

“k  g  x  3  x  ak°xx^ 


+  (Vxy}  +  FBx 


3P  .  3  , 

-  a,  - —  +  —  ( ot .  t  , 
k  9y  3  x  '  k  xy J 


(3.2) 


+  3y"  ^akayy^  +  FBy 

where  u  and  v  are  the  phasic  velocities  in  the  x  and  y  direction, 
respectively,  a  and  t  are  components  of  the  total  stress  tensor,  Fg  is 
the  body  force,  and  the  subscript  k  refers  to  the  liquid  (L)  or  vapor 
(G)  phases.  Therefore  aL  =  1-a  and  =  a. 

For  a  bubbly  flow  in  a  pipe  where  x  is  the  axial  direction  and  y  is 
transverse  or  lateral  direction  measured  from  the  wall,  the  conservation 
equations  for  the  continuous  phase,  liquid,  become 
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A  y  L  A  u  L 
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y  dir.: 


—  [(l-o)vLu.  ]  +  L-  [ ( 1-a ) Vi2] 
3x  L  3y 


.  1_  d-a)  if.  +  i_L_  [(1-a)  T  ] 
PL  9y  pL  3  x 
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p,  3y  yy 


(3.6) 


where  ML  in  the  above  equation  is  the  lift  force  on  the  bubbles  which  is 
the  only  body  force  in  the  absence  of  buoyancy. 


For  a  steady  fully  developed  flow  3/3x  =  0.  Using  the  conservation  of 
mass,  equation  (3.4),  it  can  be  shown  that  vL  =  0  and  therefore  the 
inertial  terms  in  equation  (3.6)  are  neglected.  Considering  that  the 
velocities  in  turbulent  flow  consist  of  a  mean  and  fluctuating  component 


u  =  u  +  u' 


v  =  v  +  v 


row 


The  total  stress  tensor  is  broken  down  into  viscous  or  laminar  and  tur¬ 
bulent  stresses  as  given  below. 


2 v  -  2/3  u9  -  pu'2 
3  x 


U  (|U  +  p.)  -  p7v 

3y  3  x 


,3  U  .  3  v  >  — i — i 

p  —  +  -  pu  v 

3y  3  x 


2u  ll-  -  2/3  u9  -  pv'2 


(3.7) 


_  3  u  +  3  v 
3x  3y 


In  transverse,  y  direction,  the  viscous  stresses  can  be  neglected  and 
equation  (3.6)  reduces  to 


-  (1-cO  -  (1-a)  L-  (pL  v-2)  -  pl  V'2  L-  (1-ot)  -  Ml  =  0  (3.8) 

The  conservation  of  momentum  for  the  discontinuous  phase,  vapor,  reduce1 
to  the  balance  between  the  pressure  and  lift  forces. 


(3.9) 


eliminating  3  P/3  y  from  equations  (3.8)  and  (3.9)  results  in  the 
following  equation  for  the  void  distribution  in  y  direction. 
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(3.10) 


0  1 1-3  Z  +  -i-3 
3y  ~7t 
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.,2 


In  order  to  solve  for  the  velocity  distribution  which  is  needed  for 
evaluating  the  lift  force  and  calculating  the  pressure  drop,  the  rela¬ 
tion  for  the  total  shear  stress  is  employed: 


T 


(3.11) 


The  turbulent  stress  in  the  above  equation  can  be  expressed  in  terms  of 
the  eddy  diffusivity,  non-dimensional ized  and  solved  for  the  liquid 
velocity,  if  the  eddy  diffusivity  relation  in  two-phase  flow  can  be  spe¬ 
cified.  However,  the  experimental  results  of  Refs.  (_34)  and  {25)  have 
shown  that  the  turbulent  stresses  are  strongly  affected  by  the  presence 
of  bubbles.  Therefore,  the  single-phase  relations  for  the  eddy  dif¬ 
fusivity  have  to  be  modified  to  account  for  the  turbulence  induced  by 
bubbles . 


3.1.1  Lift  Force 

The  lift  force  acts  on  a  particle  within  the  fluid  shear  layer  or  on  a 
rotating  particle  in  the  absence  of  shear.  At  high  relative  velocities, 
a  rotating  particle  may  experience  a  large  lift  force  which  is  called 
the  Magnus  effect  (36_).  This  force  acts  perpendicular  to  the  direction 
of  the  flow  and  in  the  case  of  pipe  flow  drives  the  particle  to  regions 
of  higher  velocity,  center  of  the  tube.  This  force  is  normally 
neglected  in  earth  gravity  conditions  where  the  drag  and  buoyancy  forces 
are  dominant.  At  small  relative  velocities,  the  lift  force  due  to  free 
particle  rotation  was  given  by  Rubinow  and  Keller  (37)  as 


r,  =  —  a  P  »  u 
L  8 


(3.12J 


where  fl  is  the  angular  velocity  which  for  a  freely  rotating  particle  is 
given  by 

n  =  1/2 


d  is  the  particle  diameter,  and  p  is  the  density  of  the  fluid  around  the 
particle. 

Saffman  (38),  (39)  has  shown  that  when  the  relative  velocity  between  the 
particle  and  the  liquid  is  small  and  the  fluid  Reynolds  number  is  large, 
lift  force  due  to  particle  rotation  is  less  by  an  order  of  magnitude 
than  that  due  to  shear.  The  lift  force  due  to  shear  is  given  by 


Fl  .  6.46  p1'2  (6L-UV)  2  & 1/2 


(3.13) 


where  u  is  the  fluid  viscosity. 


Lift  force  per  unit  volume  of  total  flow  become^ 
\  .  1.4  a  PL‘/2  (Vuv)  d2  (^)'/2 


(3.14) 


In  the  absence  of  gravity,  the  bubbles  will  travel  at  the  same  velocity 
as  the  local  fluid.  However,  if  the  center  of  mass  of  bubble  travels  at 
the  velocity  uL,  the  velocity  of  the  liquid  at  the  surface  of  the  bubble 
would  be 


This  means  that  the  relative  velocity  is  proportional  to  1/2  du^/dy  d. 
Assuming  that  the  Saffman  relation,  equation  (3.14),  can  be  used  for  the 
lift  force  on  the  bubbles,  the  force  per  unit  volume  becomes 


Ml  -  0.7  a  P[l/2  {^)3/2 


(3.15) 


Since  the  bubble  is  actually  surrounded  by  a  two-phase  mixture  of  vapor 
and  liquid,  an  effective  viscosity.  Ref.  (£0),  is  used  in  the  present 
analysis  to  account  for  the  existence  of  the  neighboring  bubbles. 


H,  A) 3/2 


(3.16) 


uef  f  = 


x  pl  uG  +  (1-cQ  pg  uL 
x  pl  +  (1-x)  pQ 


3.1.2  Turbulent  Eddy  Diffusivity 

Early  efforts  for  analyzing  the  velocity  distribution  in  bubbly  two- 
phase  flow  treated  the  mixture  as  a  continuous  medium.  Turbulent 
exchanges  of  momentum  and  density  were  introduced  to  evaluate  the  velo¬ 
city  and  density  (or  void)  distributions.  Levy  (4_1_)  used  this  approach 
with  the  single-phase  turbulent  mixing  length  relations  to  calculate  the 
pressure  drop  and  void  distribution.  Bankoff  (£2)  assumed  uniform  stress 
distribution  and  employed  the  single-phase  mixing  length  relation. 

Starting  with  the  original  studies  of  Serizawa  et  al.  (_35)  and  Sullivan 
et  al .  (43)  on  the  turbulence  structure  of  bubbly  two-phase  mixtures,  a 
series  of  efforts  followed  which  concentrated  on  the  effect  of  bubbles 
on  the  liquid  phase  turbulence  and  the  void  distribution  in  bubbly  flow 
(44,  45,  34) .  These  studies  have  shown  that  the  turbulence  intensity 
increases  when  bubbles  are  introduced  in  the  flow  and  in  comparison  to 


single-phase  flow  it  becomes  flatter  within  the  core.  For  moderate 
liquid  velocities,  the  turbulent  stresses  increase  with  increasing 
superficial  velocity  of  vapor,  Jq. 

In  order  to  account  for  the  effect  of  bubbles  on  the  turbulent  stresses, 
Ref.  (46)  suggests  an  eddy  diffusivity  due  to  bubble  agitation. 
Equation  (3.11)  for  the  total  shear  stress  becomes 


— - p,  u  v  -  p,  u  v 

dy  L  L 


(3.17) 


The  second  term  in  the  above  equation  is  the  single-phase  turbulent 
stress  component  and  a  single-phase  eddy  diffusivity  relation  can  be 
used  for  the  purpose  of  the  present  study.  The  third  term  on  the  right 
hand  side  is  the  enhancement  due  to  the  presence  of  bubbles.  Based  on 
the  analysis  of  potential  flow  over  a  sphere,  a  relation  for  u‘V  is 
developed  in  Ref.  (£6)  as  follows: 


u"v"  =  -  C  « 


where  C  depends  on  the  average  size  and  velocity  of  the  bubbles, 
Equation  (3.17)  will  therefore  become 


du.  du,  du. 

■  “L  -57  *  °L  e'  *  PL  c  “ -*7 


(3.18) 


For  the  single-phase  eddy  diffusivity,  Spalding's  relation,  Ref.  (47), 
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;•  =  A  (0.044)  [e°*4u+  -1  -0.4u+  -0.08u+2] 


(3.19) 


which  is  continuous  from  the  sublayer  into  the  core  is  used.  Since  this 
diffusivity  results  in  a  non-zero  velocity  gradient  at  the  center  of  the 
pipe,  Reichardt's  relation.  Ref.  (_48),  which  is  proposed  for  the  center 
regions  of  the  pipe  can  be  used  in  the  turbulent  core.  This  relation  is 


given  as 


-  —  (0.067)  [1  -  (1  -  JL)  ]  ]1  +  2  (1  -  £-)  ] 


(3.20) 


where  RQ  is  the  radius  of  the  tube. 


3.2  NON-DIMENSIONAL  EQUATIONS  AND  METHOD  OF  SOLUTION 

Substituting  the  relation  for  lift  force  from  equation  (3.16)  in 
equation  (3.10)  results  in 


-  „  ,  1/2  -1/2 

3  (1-q)  +  Ua  dv'2  +  °'7  ueff  PL  (duL^3/2  =  Q 
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v  v 
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(3.21) 


Considering  that 


T  =  Tw  d  -£-) 


equation  (18)  becomes 
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(3.22) 
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Equations  (3.21)  and  (3.22)  are  non-dimensional  ized  in  terms  of  the 
following  variables 


PL 

to  result  in 


d(l-a)  x  1-g  dv1'1'2  x  0.7  ^di /j3/2  a 
dy+  7TT2  dy+  ~+2  dy+ 


(3.23) 


(3.24) 


+ 

1  -X. 

_  _ 

1  +  C  a  +  0.044  [e°*4u+  -1  -0.4u+  -0.03u+2] 


(3.25) 


In  the  above  equation  u+  is  the  non-dimensional  mean  liquid  velocity, 
u.  \  and  the  bar  and  subscript  L  are  dropped  from  now  on.  The  above 

L  +4*2 

differential  equations  can  be  solved  simultaneously  for  a  and  u  if  v' 
variation  is  known. 

Comprehensive  measurements  of  single-phase  turbulence  structure 

including  velocity  profile,  turbulent  fluctuations,  and  Reynold  stresses 

have  been  reported  by  Laufer  (46).  Variation  of  Vv,+2  for  different 

— T3 

Reynolds  numbers  is  shown  in  Figure  3.1.  It  is  apparent  that  v'  peaks 
near  the  wall  and  then  decreases  toward  the  center  of  the  tube  and  the 
peak  value  is  a  function  of  the  Reynolds  number.  For  a  bubbly  two-phase 


flow,  the  variation  of  v'  is  similar  to  Figure  3.1,  in  that  it  peaks 
to  a  maximum  near  the  wall  and  levels  off  at  positions  aiso  close  to  the 
wall.  Refs.  (34,  35.)*  However,  measurements  of  Wang  (34)  and  Serizawa 
(35)  have  shown  that: 

1.  v'  peaks  closer  to  the  wall  in  comparison  to  single-phase  flow. 

2.  v'  is  considerably  flatter  for  two-phase  bubbly  flow. 

3.  The  maximum  value  is  a  function  of  Reynolds  number  and  quality. 

Although  the  general  trend  of  v,+2  can  be  obtained  from  the  data  of 
Refs.  (34)  and  (35),  systematic  measurement  near  the  wall  is  not 
available.  Closer  look  at  the  data  has  shown  that  detailed  variation  of 
v'  is  not  necessary  for  the  present  analysis,  and  only  the  variation 
of  maximum  v,+2  with  Reynolds  number  and  quality,  along  with  an  approxi¬ 
mate  shape  representing  the  peaking  and  flattening  of  this  parameter 
over  the  tube  cross  section,  is  sufficient. 

From  equation  (3.24)  it  can  be  seen  that  peaking  of  v,+2  will  result  in 
a  void  fraction  peak  near  the  wall.  Detailed  cross  sectional  void 
measurements  have  shown  such  a  void  peaking  for  vertical  upward  bubbly 
flow  at  earth  gravity.  Refs.  (34,  £5,  50). 

Due  to  lack  of  sufficient  data  on  the  variation  of  v,+2  in  a  bubbly  flow 
and  the  diffuculty  in  measuring  the  turbulent  fluctuations  near  the 
wall,  an  approximate  shape  for  v,+2  is  used  in  the  present  analysis. 

It  is  assumed  that  v,  +  2  rises  exponentially  near  the  wall.  It  is 
further  assumed  that  the  peak  v'  occurs  at  y/R  -  0.08.  The  following 
relations  are  used  to  specify  an  exponentially  rising  and  falling  shape 
for  the  parameter  v'+2. 


CM 


30  (JL+0.02)  -30  (i— +0.02) 

R+  R+ 


1+2  =  -J1  [-0.54  +  ^ - ]  for  <  0.08 
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R+ 


for  >  0.08 
R+ 


(3.26) 


It  should  be  emphasized  that,  although  the  relations  given  in  equation 
(3.26)  are  arbitrarily  developed  but  seem  to  be  very  specific,  the  exact 
form  of  the  equation  is  not  important.  These  equations  are  written  in 
the  above  form  to  show  an  exponentially  rising  v'+2  to  a  maximum  value 
at  y/R  =  0.08  and  an  exponentially  falling  v,+2.  Since  the  value  of 

i 

v^  is  correlated  from  our  experimental  data,  different  shapes  can  be 
used  and  v1  correlation  will  vary.  The  only  difference  may  be  that  the 
peaking  of  void  may  occur  at  a  different  radial  location. 

Equation  (3.26)  in  differential  form  becomes 


30  (JL  +  0.02)  -30  (*-  -0.08)1  2 


+  -300  (£-  -  0.08)2 

-100  (JL.  -  0.08)  e  R+ 
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X-  <  0.08 

R+ 


JL  >  0.08 
R+ 


(3.27) 


Differential  equations  (3.24),  (3.25),  and  (3.27)  can  be  solved  for  a 
and  u+  for  a  given  value  of  parameter  C. 

For  a  pipe  with  radius  R,  the  average  mass  flux  is  defined  as 
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A  =  it  R  ;  dA  =  2tt  rdr 


R 

2  /  pu  rdr 

G  = - - 

R2 

Assuming  that  liquid  and  vapor  travel  at  nearly  equal  velocities  (local 
homogeneous  condition),  the  relation  for  mass  flux  will  become 

R 

2  /Q  [«PV  +  (1-a)  PL]  u  rdr 


Non-dimensional  izing  the  above  equation  results  in 

R+ 

—  J0  [apv  +  (1-a)  PL]  u+  ( R+  -  y+)  dy+  =  y  (3.28) 
PL 


The  two-phase  Reynolds  number  is  defined  as 

Re  =  =  2R+  y  (3.29) 

ul 

The  two-phase  friction  multiplier  is  defined  as 


.  (dP) 


(3.30) 


where  the  denominator  in  the  above  equation  is  the  frictional  pressure 
gradient  if  the  total  flow  is  assumed  to  flow  as  liquid,  e.g., 


-  (lR)  =  f  JL_ 

dx  i  0  4p^R 


(3.31) 


where  f  is  the  single-phase  liquid  friction  coefficient  given  by  the 
following  relations: 


f  -  11 
0  Re 

fQ  =  1.84  x  10'4  Re0,646 


fQ  =  0.316  Re'0,26 


The  wall  shear  is  given  by 


Re  <  2300 


2300  <  Re  <  4000 


Re  >  4000 
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Therefore,  the  two-phase  friction  multiplier  becomes 


j>  2  =  — — 

0  y2  f. 


The  flow  quality  is  defined  as 


x  =  %  =  % 
X  '  w  Ga 


Wq  /  a  u  pq  dA  2pq  /0  a  u  rdr 


(3.32) 


(3.33) 


(3.34) 


(3.35) 


The  above  equation  in  non-dimensional  form  becomes 
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I o  a  u+  (R+  -  y+)  dy+ 


(3.36) 


Using  equations  (3.28)  and  (3.36)  in  (3.35)  results  in  the  following 
relation  for  flow  quality. 


x  = 


2pg  f0  a  u+  (R+  -  y+)  dy+ 


D+  2 

P,  R  Y 


(3.37) 


The  cross  sectional  average  void  fraction  is  defined  as 

R 

r  a  a  2  L  a  r  dr 
—  a  dA  J  0 
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in  non-dimensional 
R+ 

-  2  /0  «  (R 

a  =  - 
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form  this  equation  becomes 

+  +»  ,  + 

-  y  )  dy 

+  2 


(3.38) 


Equations  (3.28),  (3.37),  and  (3.38)  can  be  expressed  in  differential 
form  as  follows: 

^7  =  — 7  [a  Pv  +  (1-a)  PL]  u+  (1  -  2L-)  (3.39) 
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(3.41) 


For  a  given  value  of  the  parameters  C  and  R+,  the  differential  equations 
(3.24),  (3.25),  (3.27),  (3.39),  (3.40),  and  (3.41)  can  be  solved 

simultaneously  from  the  wall  to  R+  to  get  values  for  a,  x,  <t>Q2  and  Re 
(or  G  for  known  R).  For  given  fluid  properties  and  pipe  diameter,  there 
is  only  one  combination  of  C  and  R+  that  results  in  the  desired  G  and  x 
combination.  Therefore,  for  every  G  (or  Re)  and  x  combination,  the  l 

solution  for  void  fraction  and  t^o-phase  friction  multiplier  can  be 
obtained. 


3.2.1  Boundary  Conditions 


In  order  to  solve  the  above  differential  equations,  the  boundary  con¬ 
ditions  at  the  wall  should  be  known.  Since  the  non-dimensional  tur- 
bulent  fluctuation,  v'  is  zero  at  the  wall,  equation  (3.24)  is  not 
defined  at  the  wall.  The  integration  starts  from  a  point  very  close  to 
the  wall  where  all  liquid  conditions  can  be  assumed  and  the  state 
variables  can  be  defined.  y+  =  0.1  is  within  the  laminar  sublayer  and 
sufficiently  close  to  the  wall  where  a  zero  void  assumption  can  be  made. 
Therefore,  at  y+  =  0.1 
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u  =  0.1 


Y  =  2_CIMl!.iO (3.42) 

R+  ?  3R  R 


x  =  0 


a  =  0 


3.2.2  Method  of  Solution 

A  computer  program,  ZGBUB,  was  developed  to  solve  the  six  differential 
equations,  (3.24),  (3.25),  (3.27),  (3.39),  (3.40),  and  (3.41)  from  the 
boundary  condition  given  above  up  to  a  specified  R+.  This  program  uses 
the  SFODE  integration  algorithm  which  solves  a  system  of  ordinary  dif¬ 
ferential  equations  using  Adams  Predictor-Corrector  method.  ZGBUB  was 
written  in  FORTRAN,  and  the  calculations  were  performed  on  an  IBM-PC/XT 
(with  MS-DOS  operating  system)  and  a  MicroVax-II  (with  VMS  operating 
system).  All  the  routines  were  written  in  double  precision  to  reduce 
the  round-off  errors. 

In  order  to  solve  the  differential  equations  for  a,  x,  <f>  2t  and  G  the 
values  of  v'm  in  equation  (3.26)  and  parameter  C  in  equation  (3.25) 
should  be  known.  The  experimental  data  for  silicone-water  and  Dowtherm- 
water  experiments  described  in  Section  2  was  used  to  obtain  the  values 
of  C  and  v'm.  As  will  be  discussed  later,  only  a  correlation  for  v'  in 
terms  of  Re  and  x  is  needed  to  close  the  problem.  For  a  given  set  of 
fluid  properties  and  pipe  diameter,  there  are  only  unique  values  of  C 
and  R+  which  can  result  in  a  desired  G  and  x  combination.  Therefore, 
iteration  on  different  values  of  C  and  R+  can  be  performed  to  give  the 
needed  G  and  x  which  will  also  provide  $Q2  and  a  values. 


3.3  RESULTS  AND  DISCUSSION 

The  empirical  parameters  in  the  present  analysis  are  v1  which  is  the 
maximum  value  for  the  turbulent  intensity  (v,+2)  and  C  which  is  the 
coefficient  for  the  eddy  diffusivity  due  to  the  presence  of  bubbles.  As 
mentioned  earlier,  the  data  from  Refs.  (34) ,  (44) ,  (45)  and  (46_)  shov,* 
that  the  turbulent  intensity  peaks  to  a  maximum  value  near  the  wall  (y/R 
between  0.05  to  0.1)  and  unlike  the  single-phase  flow,  it  decreases 
rapidly  and  flattens  out  for  y/R  >  0.1.  In  addition,  the  value  of  v'  is 
a  function  of  liquid  Reynolds  number  and  the  vapor  phase  content.  The 
parameter  C  depends  on  the  average  bubble  size  and  velocity. 

The  experimental  pregram  and  results  for  silicon-water  and  Dowtherm- 
water  tests  were  discussed  in  the  previous  sections.  The  reduced  data 
consisted  of  cross  sectional  average  void  fraction  and  two-phase  fric¬ 
tion  multiplier  for  the  given  mass  flux  (G)  and  quality  (x).  As  men¬ 
tioned  earlier,  due  to  high  phase  velocities,  photographic  methods  for 
visual  determination  of  flow  regimes  were  not  possible  in  the  present 
experiments.  However,  the  differential  pressure  measurements  and  the 
Auburn  meter  measurement  gave  an  indication  of  departure  from  bubbly 
flow.  Generally,  this  occured  around  void  fraction  of  0.3.  Therefore, 
the  tests  where  we  had  a  bubbly  flow  simulation  could  be  identified  and 
used  in  the  present  analysis. 

Silicon-water  and  Dowtherm-water  properties  were  used  in  the  ZGBUB 
program  to  find  the  measured  two-phase  friction  multiplier  for  the  given 
combination  of  mass  flux  and  flow  quality.  Only  the  bubbly  flow  tests 
were  used  for  this  purpose  and  the  test  conditions  are  given  in  Table 
3.1.  Only  one  combination  of  R+,  v'  and  C  could  be  found  to  match  the 
data.  The  predicted  values  of  v'  and  C  are  given  in  Table  3.2. 

A  computer  program  which  obtains  a  weighted  least  square  fit  of  a  user 
specified  function  to  data  values  by  means  of  stepwise  Gauss-Newton 
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Table  3.1 


TEST  CONDITIONS  FOR 

BUBBLY  FLOW  SIMULATION  TESTS 
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TEST  # 

Volume  Flow 
(6PM) 

Silicon  Water 

Total  Mass  Flux 
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54 
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£ 
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Table  3.2 


VALUES  OF  v'  AND  C 

FOR  BUBBLY 

FLOW  SIMULATION  TESTS 

Total  Mass  Flux 

Qua! i ty 

v' 

C 

(Lbm/ft2/s) 

767.1 

0.083 

5460 

-18.40 

767.0 

0.295 

19000 

-3.39 

767.0 

0.259 

14100 

-5.79 

766.6 

0.189 

7900 

-6.15 

766.6 

0.153 

6700 

-8.60 

767.1 

0.083 

5310 

-  14.00 

767.0 

0.295 

18000 

-2.70 

852.6 

0.127 

4775 

-5.70 

852.1 

0.190 

5730 

-3.00 

852.5 

0.318 

12950 

-1.10 

853.0 

0.095 

4629 

-10.20 

852.5 

0.159 

5158 

-3.90 

1023.6 

0.075 

3184 

-3.45 

1023.9 

0.260 

4775 

1.00 

1024.0 

0.155 

3811 

-1.80 

1122.6 

0.109 

2994 

-0.50 

1123.0 

0.181 

3280 

1.55 

1123.4 

0.278 

4224 

2.10 

899.8 

0.072 

1776 

12.00 

898.7 

0.139 

2027 

9.55 

897.6 

0.205 

2305 

7.10 

796.9 

0.086 

2107 

6.30 

795.7 

0.161 

2412 

3.50 

794.6 

0.236 

2940 

1.00 

449.8 

0.157 

4730 

-3.80 
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iterations  was  used  to  correlate  v1  and  C.  This  nonlinear  least 
squares  regression  program,  BMD07R,  is  part  of  the  SIGSTAT  statistical 
package.  Ref.  (52.),  and  can  have  as  many  as  10000/n  variables,  where  n 
is  the  number  of  observations.  For  our  purposes  there  are  two  indepen¬ 
dent  variables  namely,  Reynolds  number  defined  as 

Re  =GD 

and  quality  x. 

The  predicted  values  for  v'm  and  C,  given  in  Table  3.2,  were  used  to 
find  the  functions  of  independent  variables  which  represent  the 
variations  of  v‘m  and  C.  These  functions  were  then  used  in  BMD07R  to 
find  the  parameters  which  would  fit  the  predicted  values.  The  best  fit 
for  v1  values  is  given  below. 

-5.0[x”0*387  -7 . 7 ( Re  lCf3)’0*499]  -1.038 

v'm  =  34.71  e  +  345.95(10"5Re) 

(3.43) 

Comparison  of  the  predicted  values  of  v'm  from  Table  3.2  and  those 
obtained  from  the  above  relation  is  shown  in  Figure  3.2.  The  deviation 
is  within  ±20%.  A  fit  for  the  parameter  C  in  terms  of  the  Reynolds 
number  and  quality  was  also  found.  It  should  be  noted  that  since  C 
depends  on  the  average  bubble  size  and  velocity,  this  type  of  fit  will 
not  be  useful  to  define  its  variation  and  a  third  independent  parameter 
may  be  needed.  This  parameter  should  depend  on  the  fluid  properties 
including  interfacial  tension.  However,  data  from  a  range  of  fluids  with 
varying  properties  is  needed  to  establish  such  a  parameter. 

As  mentioned  earlier,  for  a  given  combination  of  G  and  x  there  is  a 
unique  solution  which  can  be  obtained  with  specific  values  of  C  and  R+ 
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Therefore,  functional  variation  of  the  parameter  C  is  not  necessary  for 
obtaining  the  predictions  for  the  two-phase  friction  multiplier  and  void 
fraction.  However,  knowledge  of  such  a  correlation  will  eliminate  the 
need  for  iteration  on  C  and  R+  values  in  order  to  match  the  G  and  x  com¬ 
bi  nation. 

Obviously,  the  two-phase  friction  multiplier  predicted  by  the  present 
method  would  fit  the  silicon-water  and  Dowtherm-water  data  fairly  well 
since  it  was  based  on  those  results.  A  comparison  of  the  predicted  and 
measured  two-phase  multipliers  is  shown  in  Figure  3.3.  It  should  be  noted 
that  the  data  and  predictions  for  void  fraction  and  friction  pressure 
drop  are  very  close  to  the  results  of  the  Homogeneous  Equilibrium  Model 
(HEM).  The  two-phase  friction  multiplier  and  void-quality  relation  from 
HEM  are: 


=  1  + 


>L  '  PG 


(3.44) 


(3.45) 


+  ^  (1-x) 


The  present  model  assumes  a  local  homogeneous  condition  but  allows  for 
void  distribution  within  the  cross  section.  This  results  in  void  peaking 
near  the  wall  but  nearly  a  flat  void  distribution  a  small  distance  away 
from  the  wall.  The  reason  for  homogeneous  behavior  of  the  data  is 
believed  to  be  the  equal  density  of  fluids  rather  than  an  inherent  con¬ 
dition  in  the  absence  of  buoyancy.  If  the  model  developed  based  on  the 
above  data  is  mechanistically  based  and  represents  the  physical  phenome¬ 
non,  it  should  account  for  different  fluid  conditions  and  will  not 
necessarily  predict  close  to  HEM  at  all  times.  Actually,  it  will  be 
shown  later  that  the  proposed  model  predicts  considerably  different 
pressure  drop  behavior  than  HEM  for  gas-liquid  flow  at  zero  gravity. 
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However,  the  predicted  cross  sectional  average  void  fractions  are  essen¬ 
tially  the  same  as  HEM.  Apparently,  the  void  peaking  near  the  wall  and 
the  absence  of  phase  slip  results  in  much  larger  two-phase  pressure 
drop. 

In  order  to  evaluate  the  predictive  capability  of  this  model,  it  should 
be  compared  with  the  true  zero  gravity  adiabatic  gas-liquid  data.  As 
discussed  earlier,  there  is  practically  no  zero  gravity  bubbly  flow  data 
which  can  be  used  for  our  purpose.  Results  of  a  series  of  adiabatic  two- 
phase  flow  tests  aboard  a  KC-135  aircraft  flying  in  a  parabolic  trajec¬ 
tory  have  been  reported  in  Ref.  (52).  These  tests  were  mostly  performed 
in  slug  and  annular  flow  regimes,  and  in  addition  to  the  flow  regime 
transitions,  have  qualitatively  shown  that  in  comparison  to  earth  gra¬ 
vity  conditions,  the  pressure  drop  significantly  increases.  As 
discussed  earlier,  the  majority  of  two-phase  flow  studies  at  zero  gra¬ 
vity  are  presently  in  progress  and  the  results  should  be  available  by 
early  1988.  Among  these,  only  the  results  of  a  series  of  tests  by 
Battelle  Laboratories,  Ref.  (53),  have  been  made  available  to  us.  These 
tests  consisted  of  a  set  of  convective  boiling  tests  with  water  and  a 
set  of  adiabatic  nitrogen/water  tests  at  atmospheric  conditions. 
Unfortunately,  these  tests  did  not  cover  the  whole  range  of  flow  regimes 
and  only  two  sets  of  conditions  may  have  been  in  the  bubble  regime. 
Close  examination  of  these  tests  showed  that  the  flow  was  laminar  and 
the  present  model  can  not  be  used  to  predict  the  results.  However,  it 
was  generally  concluded  that,  (27),  the  pressure  drop  was  higher  than 
expected  under  earth  gravity  conditions.  The  air/water  tests  performed 
by  the  University  of  Houston,  Ref.  (_54 ) ,  aboard  the  NASA  learjet  covered 
the  entire  range  of  flow  regimes.  These  tests  were  mainly  aimed  at  flow 
regime  identification  and  the  pressure  drop  results  have  not  been 
published  yet. 

Since  the  present  model  does  not  have  a  closed  form  solution  and  contains 
an  iterative  approach,  solutions  can  only  be  developed  for  specific  fluid 


i 


conditions  and  test  section  geometry.  In  order  to  provide  predictions 
which  can  be  used  in  design  of  spacecraft  thermal  management  systems,  a 
series  of  results  for  selected  fluids  are  presented  here.  In  addition, 
the  predictions  for  the  test  conditions  of  the  zero  gravity  experimental 
efforts  in  progress  are  presented  and  compared  to  earth  gravity  models. 
These  predictions  can  be  compared  with  the  test  results  when  published. 

A  majority  of  the  separated  flow  models  for  two-phase  friction  pressure 
drop  under  earth  gravity  conditions  are  empirical  and  only  a  few  have 
physical  basis.  These  models  do  not  differentiate  between  different  flow 
regimes,  boiling  or  adiabatic  conditions,  and  the  test  section  orien¬ 
tation.  The  most  widely  used  model  is  the  Lockhart-Martinel 1 i -Nel son 
correlations.  The  Lockhart-Martinel 1 i  model  was  developed  for  adiabatic 
flow  of  low  pressure  air-water  mixtures  in  horizontal  tubes.  The  two- 
phase  friction  multiplier  was  correlated  as  a  function  of  the  Mart i nel 1 i 
parameter  defined  as 
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The  general  form  of  the  two-phase  friction  multiplier  (based  on  the 
single-phase  component  of  pressure  drop  calculated  for  liquid  flowing 
alone  in  the  channel )  is 


=  1  + 


(3.46) 


where  the  values  of  parameter  C  are  given  for  Laminar/turbulent  com¬ 
binations  of  gas  and  liquid  phases.  For  the  case  of  turbulent  flow  in 
both  gas  and  liquid,  C=20  and  Martinelli  parameter  can  be  approximated  by 
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Although,  this  correlation  was  developed  for  air-water  flow,  the  results 
were  presented  in  general  form  to  be  used  for  steam-water  mixtures. 

This  model  was  extended  by  Martinel  1  i -Nel son  (_66)  to  forced  convective 
boiling  of  water  as  follows.  It  was  assumed  that  the  flow  of  gas  and 
liquid  will  be  turbulent.  It  was  further  assumed  that  the  air-water 
results  of  Lockhart -Marti nel 1 i  can  be  used  for  steam-water  at 
atmospheric  conditions.  A  relation  for  <j>L2  at  critical  pressure  was 
developed  by  assuming  equal  densities  and  viscosities  for  both  phases. 
Two-phase  multipliers  for  intermediate  pressures  were  then  interpolated 
from  atmospheric  and  critical  values.  For  steam-water  mixtures  at  hiqh 
pressures,  two-phase  multipliers  are  generally  enveloped  between 
Martinelli-Nelson  and  HEM  predictions  and  approach  HEM  as  mass  flow 
rates  are  increased.  Unlike  the  high  pressure  conditions,  HEM  could 
result  in  higher  two-phase  multipliers  at  low  pressures.  The  predictions 
by  Lockhart -Marti nel 1 i  and  HEM  for  air-water  mixtures  at  atmospheric 
conditions  are  given  in  Figure  3.3.  Other  correlations  include  Thon, 
(55) ,  which  is  applicable  to  steam-water,  and  Baroczy,  ( _56 ) ,  which  is 
applicable  to  a  wide  range  of  boiling  liquids. 

For  steam-water  mixtures  at  1000  psi,  predictions  by  the  present  model 
for  G=358  Lbm/ft2-s  (1750  Kgm/m2-s)  are  compared  to  predictions  by  HEM, 
Thom,  and  Mart inel 1 i -Nel son  along  with  the  earth  gravity  data  in  Figure 
3.4.  The  effect  of  mass  flux  on  the  zero  gravity  model  in  this  range  is 
small.  In  general  it  was  found  that,  unlike  the  earth  gravity  con¬ 
ditions  where  the  two-phase  friction  multiplier  decreases  with 
increasing  mass  flux,  at  zero  gravity  the  multiplier  may  increase  or 
decrease  with  mass  flux  depending  on  the  fluid  condition.  It  is 
believed  that  decreasing  with  increasing  mass  flux  at  1  g  conditions 
is  due  to  phase  slip  which  does  not  exist  at  zero  gravity. 


When  the  mass  flux  increases,  the  slip  between  the  phases  is  reduced 
which  causes  4>0  to  decrease  and  approach  the  homogeneous  value.  Under 
zero  gravity  conditions,  there  is  no  slip  between  the  phases  and  such  a 
trend  is  not  expected  to  occur.  The  two-phase  friction  multiplier 
variation  with  mass  flux  is  dictated  by  the  relative  properties  of  the 
phases. 

The  zero  gravity  tests  performed  for  air/water  mixtures  in  a  0.5  inch 
tube  at  atmospheric  conditions,  Ref.  (54.),  are  for  mass  flow  rates  of  96 
and  186  Lbm/ft2/s.  The  predictions  by  the  present  model,  HEM,  and 
Lockhart -Martinell i  for  these  conditions  are  shown  in  Figure  3.5.  These 
results  can  be  compared  with  data  when  published. 

Ammonia  requires  considerably  lower  pumping  power  for  equal  heat  removal 

in  comparison  to  other  refrigerants  at  50-80°F  which  is  the  range  of 

operating  condition  for  the  two-phase  thermal  management  loops. 

Predictions  of  the  two-phase  friction  multiplier  for  ammonia  in  a  1  inch 

tube  at  different  pressures  and  mass  flow  rates  are  shown  in  Figures  3.6 

and  3.7.  Unlike  the  air/water  flow,  increasing  mass  flow  rates  results 
2 

in  increasing  4>0  ,  Increasing  system  pressure  results  in  decreasing 
two-phase  friction  multiplier. 

As  mentioned  earlier,  the  void  fractions  predicted  by  the  present  model 
are  slightly  below  the  HEM  predictions  but  can  generally  be  approximated 
by  the  HEM  relation,  equation  (3.45). 

3.4  BUBBLY  FLOW  WITH  NUCLEATION  AT  THE  WALL 

In  a  two-phase  thermal  transport  loop,  large  amounts  of  heat  are  removed 
from  the  heating  surface  due  to  the  latent  heat  of  vaporization  of 
liquid.  Normally,  single  phase  liquid  enters  the  heating  section  and 
nucleation  takes  place  at  the  wall.  In  comparison  to  a  single  phase 
loop,  this  process  will  reduce  the  pumping  power  and  increase  the  heat 
transfer  coefficients  by  up  to  two  orders  of  magnitude. 
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Figure  3.5  -  Two-Phase  Friction  Multiplier  vs.  Quality 
for  Air/Water  at  14.7  psi  in  a  0.5  inch  Diameter  Tub 


100  Lbm/ft2/s 
1 .  inch 


As  mentioned  earlier,  the  empirical  models  developed  for  two-phase  fric¬ 
tion  multiplier  and  void-quality  relation  at  earth  gravity  do  not  dif¬ 
ferentiate  between  adiabatic  and  diabatic  conditions  or  different  flow 
regimes.  In  order  to  see  the  effect  of  wall  voidage  on  the  pressure 
drop,  the  model  developed  for  adiabatic  conditions  at  zero  gravity  can 

be  used  and  the  effect  of  wall  voidage  can  be  incorporated  through  the 

boundary  condition.  This  means  that  for  a  given  wall  heat  flux  and  sur¬ 
face  conditions,  the  void  fraction  at  the  wall  should  be  evaluated.  It 

should  be  noted  that  the  adiabatic  model  did  not  include  any  mass  or 

momentum  transfer  across  the  interface.  This  would  obviously  be  a  major 
assumption  for  diabatic  two-phase  flow  of  a  single  component  fluid. 

It  is  known  that  the  number  and  size  of  the  bubbles  at  the  wall  during 
the  nucleate  boiling  process  depend  on  a  variety  of  parameters  including 
the  wall  heat  flux,  fluid  properties,  surface  condition,  and  mass  flow 
rate  (for  convective  boiling  process).  Therefore  determination  of  the 
void  fraction  at  the  wall  in  terms  of  the  above  factors  can  not  be 
easily  accomplished.  The  approach  adopted  here  is  to  estimate  the  void 
fraction  for  a  specific  condition,  i.e.  critical  heat  flux,  and  to 
correlate  the  void  fraction  at  the  wall  in  terms  of  the  known  parame¬ 
ters. 

It  is  known  that  the  number  of  active  nucleation  sites  is  proportional 
to  the  heat  flux  according  to  the  following  relation 

.  n 

N  =  D  (q) 


Empirical  information  on  a  wide  variety  of  liquids  and  surfaces  has 
shown  that  the  range  of  variation  of  the  exponent  n  is  between  2  and  3. 
This  does  not  imply  that  the  heat  transfer  is  independent  of  surface 
finish  since  the  coefficient  D  is  a  function  of  the  surface  conditions. 
For  a  given  surface  with  known  cavity  distribution  it  can  be  concluded 
that  the  void  at  the  wall,  aw>  -js  related  to  heat  flux  by  the  above 
relation 

3-34 


v -v-*  .• -/ . 


-V 


1 


2 


v^v.v.w  v  w.v  w.vw  vy 


°w  =  E  ^ 


Therefore,  if  the  critical  heat  flux  and  the  corresponding  void  fraction 
at  the  wall  can  be  determined,  they  can  be  used  to  find  the  void  for  any 
other  heat  flux  by 


“w  “  ^CHF 


(-ri-) 


2 


qCHF 


where  q^p  and  («w)(;np  are  the  critical  heat  flux  and  wall  void  fraction 
respectively. 

It  is  generally  agreed  that  the  mechanism  leading  to  critical  heat  flux 
(CHF)  depends  on  the  quality.  At  high  qualities  when  the  flow  is  annular 
CHF  is  due  to  dryout  of  the  liquid  film  and  should  be  modeled  by 
balancing  the  evaporation  and  entrainment  with  droplet  deposition.  At 
low  qualities,  typical  of  bubbly  flow  conditions,  CHF  is  a  local  pheno¬ 
menon  and  is  caused  by  build  up  of  a  bubbly  layer  at  the  wall.  Ref. 
(57).  As  bubble  nucleation  increases,  the  lift  and  turbulent  forces  are 
insufficient  to  transport  the  bubbles  to  the  regions  of  higher  velocity. 
Therefore,  the  bubble  layer  will  become  thicker  and  at  CHF  when  it  has 
reached  the  maximum  thickness  the  enthalpy  transport  with  the  core  is 
limited.  With  this  mechanism  leading  to  the  CHF  condition  the  void  frac¬ 
tion  at  the  wall  can  be  readily  estimated. 

If  the  average  size  of  the  bubbles  attached  to  the  wall  is  known,  the 
local  void  fraction  can  be  determined  from  the  following  relation 

a  =  (Number  of  Bubbles)  (Volume  of  a  Bubble) 

Total  Volume 

If  the  bubbles  on  the  wall  are  separated  by  a  distance  of  S,  the  number 
of  bubbles  along  the  tube  is  given  by  Ph/s,  where  Ph  is  the  perimeter  of 
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the  tube.  The  above  relation  for  spherical  bubbles  with  average  radius 
of  rc  would  therefore  become 


°w  = 


(it  D/S)  (4/3  it  rc3) 
n  D  S  r 

c 


(3.48) 


The  average  size,  r  ,  does  not  affect  the  present  analysis  but  if  the 
radius  at  the  bubble  departure  is  r^,  it  can  be  assumed  that  the  average 
radius  of  the  attached  bubbles  is  r^/2. 

It  is  postulated  that  CHF  occurs  for  the  maximum  void  fraction  at  which 
individual  bubbles  can  be  maintained.  Therefore,  aw  can  be  estimated 
from  the  condition  that  the  bubbles  would  be  just  touching  one  another, 
i.e.  S  =  2  rc.  The  local  wall  void  fraction  at  CHF  from  the  above 
equation  becomes  («w)qHp  =  0.52.  According  to  Ref.  (57) ,  photographic 
studies  of  the  bubbly  layer  have  shown  that  the  shape  of  the  bubbles  can 
be  approximated  by  an  ellipsoid  with  ratio  of  long  to  short  axis  of  3:1. 
With  such  bubbles  just  touching  the  void  fraction  at  CHF  becomes 
( aw)cHF  =  0,82  and  the  loca1  wa^  void  fraction  will  be 


%  =  0.82  (-S-9-)  (3.49) 

qCHF 


The  problem  of  establishing  the  critical  heat  flux  under  pool  boiling 
condition  at  zero  gravity  has  been  studied  extensively  and  reported  among 
others  in  Refs.  (58)  (59)  (60)  (61).  It  is  generally  believed  that  the 
critical  heat  flux  decreases  with  gravity  as  g*^.  This  corresponds  to 
the  earth  gravity  theories  of  critical  heat  flux  and  has  been  confirmed 
for  gravity  range  of  0.01  to  1.0  g  (58)  ( _59) .  However  drop  tower  tests 
of  Sherley  (62)  and  tests  with  magnetic  fluids,  (_60) ,  have  shown  that  such 
a  dependence  on  gravity  breaks  down  for  g  <  0.1  gQ  ancj  the  critical  heat 
flux  approaches  a  value  corresponding  to  0.5  to  0.7  of  the  earth  gravity 
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CHF.  Based  on  the  existing  data,  it  can  be  concluded  that  no  satisfac¬ 
tory  conclusion  regarding  the  pool  boiling  CHF  at  zero  gravity  can  be 
drawn.  There  is  practically  no  information  available  on  the  critical 
heat  flux  under  two-phase  convective  boiling  conditions.  The  study 
reported  in  Ref.  (63)  has  shown  that  the  effect  of  buoyancy  on  critical 
heat  flux  at  earth  gravity  decreases  with  increasing  pressure,  inlet 
velocity,  and  subcooling.  It  can  be  assumed  that  the  effect  of  gravity 
on  decreasing  the  critical  heat  flux  would  be  much  smaller  than  the  pool 
boiling  case,  since  liquid  drag  will  cause  the  bubbles  to  depart  in  con¬ 
vective  boiling  conditions.  Recent  modeling  efforts  for  earth  gravity 
forced  convective  critical  heat  flux  are  reviewed  in  Ref.  (64)  and 
results  of  several  experimental  studies  are  shown  in  Ref.  (65). 

In  the  absence  of  a  correlation  for  the  critical  heat  flux,  the  effe’ct 
of  wall  voidage  can  be  qualitatively  shown  for  different  heat  fluxes. 
Using  equation  (3.49)  for  q/qCHF  of  0.05  and  0.15  the  void  fraction  at 
the  wall  would  be  0.002  and  0.018,  respectively.  The  model  developed  for 
the  adiabatic  conditions  was  used  for  refrigerant  11  with  wall  void 
fraction  of  zero,  0.002,  and  0.018.  The  results  are  shown  in  Figure  3.8 
and  it  shows  less  than  15%  increase  in  the  two-phase  friction  multiplier 
for  q/qCHF  of  0.15  and  negligible  Increase  for  q/qCHF  of  0.05.  Since  the 
proposed  two-phase  loops  for  the  spacecraft  cooling  systems  are  expected 
to  operate  at  low  power  levels  (~1  W/cm2),  the  effect  of  wall  voidage 
would  be  small  unless  the  critical  heat  fluxes  for  the  operating  con¬ 
ditions  are  found  to  be  smaller  than  5  W/cm2. 
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Section  4 


ANALYSIS  OF  ANNULAR  FLOW 


It  is  generally  accepted  that  under  annular  flow  conditions,  a  signifi¬ 
cant  fraction  of  liquid  may  be  entrained  as  droplets  in  the  gas  core. 
There  is  continuous  interchange  of  mass  and  momentum  between  the  liquid 
film  and  the  gas  core.  The  predictive  capability  of  a  phenomenologi¬ 
cally  based  annular  flow  model  rests  on  the  relations  used  to  define  the 
interfacial  friction  factor  and  the  entrainment  fraction. 

In  this  section,  a  one-dimensional  steady  state  annular  flow  is  pre¬ 
sented.  The  model  uses  the  triangular  relationship  proposed  by  Hewitt 
(31) .  By  using  the  experimental  results  discussed  in  Section  2,  empiri¬ 
cal  relations  for  the  interfacial  friction  factors  and  entrainment  frac¬ 
tion  are  obtained.  The  resulting  model  provides  a  method  for  evaluating 
the  two-phase  friction  multiplier  and  void  fraction  under  annular  flow 
conditions  in  the  absence  of  buoyancy. 

4.1  ONE-DIMENSIONAL  STEADY-STATE  ANNULAR  FLOW  MODEL 

Consider  a  horizontal  steady  fully  developed  annular  flow  as  shown  in 
Figure  4.1.  Assuming  that  the  core  can  be  considered  as  a  homogeneous 
mixture  of  gas  and  entrained  liquid,  the  mean  core  density  is  given  by 

w 

°c  “JT-G  ■  »G  <>  <4-l> 

G  x 

where  e  is  the  entrainment  fraction  and  w  is  the  mass  flow  rate  in  the 
gas  core 

wc  =  wr  +  e  w. 


(4.2) 


It  is  also  assumed  that  the  ratio  of  the  mass  flux  in  the  core  to  mass 
flux  within  the  liquid  film  is  constant,  i.e. 


PG  qG  *  e  PL  \ 

- i -  (4.3) 

(1-e)  pl  qL 


where  ac  is  the  ratio  of  the  core  cross-sectional  area  to  the  total 
cross  sectional  area  of  the  tube 


“c=<l-^)2 


(4.4) 


Equation  4.3  can  be  solved  for  the  entrainment  fraction  which  results  in 
the  following  equation 


e  =  [ 


1  -  a. 


K  a. 

]  C  - £- 


PG  qG 


1  +  ( K-l )  ac  1  -  ac  pL  qL 


(4.5) 


The  mean  core  velocity  is  given  by: 


U.  =  — —  G  (x  +  (l-x)e) 


(4.6) 


The  relation  between  the  actual  void  fraction  a  (when  entrainment  is 
considered)  and  a  can  be  expressed  as: 


» 


The  definition  of  the  interfacial  friction  factor  can  be  written  as: 


2  c  c 


(4.8) 


A  force  balance  on  the  gas  core  in  the  z  direction  gives  (see  Figure 
4.1): 


t,  =  _Il  {  I?  +  (l£)2  L.  (G2  (^)2  £-) 
2  3z  r.  3Z  r, 


(4.9) 


i  9 


Since  fully  developed  flow  is  considered,  (4.9)  is  simplified  to: 

T  H  if. 

Ti  =  '  2  3 1 


(4.10) 


Note  that 
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(4.11) 


which  yields  from  equation  (4.10) 


(4.12) 


The  relation  for  the  shear  stress  distribution  in  the  liquid  film  given 
by  Hewitt  (31)  is: 


(r)  =  t  (II)  .ilf  (II-— ) 

1  r  2  dz  r 


(4.13) 


Note  that 


(4.14) 


(4.15) 


Equation  (4.13)  can  be  written  as: 


T  (y)  =  x. 
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The  shear  stress  in  the  liquid  film  is  assumed  to  be  given  by  the 
relation: 
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(4.17] 


where  e  is  an  assumed  given  turbulent  eddy  diffusivity  correlation. 
With  the  above  relation,  the  film  thickness  and  thus  the  two-phase 
multiplier  and  void  fraction  can  be  determined  as  follows. 

Integrating  (4.17)  from  y  =  0  to  y  results  in  u(y)  which  can  be  used  to 
calculate  the  total  liquid  film  mass  flow  rate  by: 


m*  =  /  *(D-2y)  p£  u(y)  dy 
o  * 


(4.18) 


For  each  guessed  value  of  6,  the  calculated  m^  can  be  compared  with  the 
actual  liquid  film  mass  flow  rate: 


m*  =  *  D2G  (l-x)(l-e) 


(4.19] 


When  the  two  values  for  m  are  equal,  the  film  thickness  is  obtained. 
The  Mart i nel 1 i ' s  parameter  Xu  is  defined  in  2.5.3  and  is  writen  as 
,pi  1/2  ,  „ 


i 


The  two-phase  multiplier  <j>^  defined  in  2.5.2  is  given  by 


(dP) 

vdz;F 


f  "  © 
dzf 


4.2  N0N-DIMENS10NALIZATI0N  OF  EQUATIONS 

The  interfacial  friction  factor  in  (4.6)  can  be  expressed  as 


Ti  =  I  fipcuc2 
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where 


a  =  — 
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Similarly,  dividing  (4.16)  by  _  gives 
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where 


Re*  =  GD 


From  (4.18)  we  have 


dm. 


—  =  w  (D-2y)  P*  u(y) 


Defining 


M  2 


4m  o 


ttD2G 


and  using  (4.27)  and  (4.29)  yields 


dM 

dn 
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where 


Using  (4.19)  and  (4.29)  results  in 
M  =  (1-x)  ( 1-e) 


(4.28) 


(4.29) 


(4.30) 


(4.31) 
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The  two-phase  multiplier  in  (4.21)  can  be  related  to  the  nondimensional 
quantities  as  follows. 


/dPv  =  £[*  G2  (1-x)2 
dz  f  D  p. 


Using  (4.21),  equation  (4.24)  can  be  written  as 


(4.32) 
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(4.33) 


4.3  PROCEDURE  FOR  CALCULATING  THE  FILM  THICKNESS,  TWO-PHASE  MULTIPLIER, 
AND  VOID  FRACTION 

From  the  above  results  and  the  given  correlations  for  the  interfacial 
friction,  f^ ,  and  the  turbulent  eddy  diffusivity,  e,  of  the  liquid  film, 
the  film  thickness  can  be  obtained  by  iteration  as  described  below. 

(1)  Given  flow  conditions  G,  x  and  fluid  properties. 

(2)  Choose  values  for  dimensionless  liquid  film  thickness,  a. 

(3)  Use  the  correlation  for  the  interfacial  friction  factor,  f . ,  and 
(4.22)  to  calculate  t.*. 

(6)  Use  (4.4)  and  (4.26)  to  calculate  (— )* 

dz 

(5)  Use  the  relation  for  eddy  diffusivity  and  (4.25)  to  integrate 
(4.28)  from  n  =  0  to  n  =  1  to  obtain  u*  for  different  n  values. 


I 
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(6)  Use  the  result  of  step  5  and  (4.30)  to  calculate  _  and  integrate 

dn 

« 

from  n  =  0  to  to  n  =  1  to  obtain  the  liquid  film  mass  flow  rate,  M. 

(7)  Compare  the  result  of  step  6  with  equation  (4.31).  If  the  values 
are  equal  then  the  chosen  value  of  a  is  correct;  otherwise,  a  new 
value  of  a  is  chosen  and  the  same  procedure  is  repeated. 

(8)  With  a  correct  value  of  a,  the  Martinelli  parameter,  Xtt,  void 

fraction,  a^t  and  the  two-phase  friction  multiplier,  can  be 

obtained  from  equations  (4.20),  (4.7),  and  (4.33). 


4.4  INTERFACIAL  FRICTION,  ENTRAINMENT  RATE,  AND  EDDY  DIFFUSI VITY 
CORRELATIONS 

As  mentioned  earlier,  the  accuracy  of  the  physically  based  annular  flow 
models  rests  on  the  correlations  used  for  the  interfacial  friction  and 
entrainment  rate.  Wallis,  Ref.  (.68),  suggests  the  following  relation 
for  the  interfacial  friction  factor: 


fi  ■  U  <l  +  300  D1 


where  f$  is  the  single-phase  friction  factor  for  the  gas  flowing  alone 
in  the  tube  with  the  total  core  flow  rate.  Wallis  uses  a  value  of  0.005 
for  f  .  Whalley  and  Hewitt,  Ref.  (J59),  use  the  following  relation  for 
the  interfacial  friction  factor: 


fi  *  f 5  [  1  ♦  <!r> 


(4.34) 


Different  models  for  entrainment  rate  have  been  proposed  but  no  satis¬ 
factory  correlation  is  presently  available. 
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For  calculation  of  the  velocity  profile  in  the  liquid  film.  Ref.  (31) 
uses  Deissler's  relation,  while  Levy,  (.70}  uses  the  mixing  length 
approach.  The  model  described  in  Section  4.3  was  used  to  predict  the 
two-phase  friction  multiplier  and  cross-sectional  average  void  fraction 
for  the  conditions  tested  with  L-45/water.  Different  relations  for  the 
interfacial  friction,  entrainment  fraction,  and  eddy  diffusivity  were 
tested,  and  the  following  relations  were  found  to  result  in  best 
agreement. 

Interfacial  Friction  Factor: 

An  exponential  variation  with  the  film  thickness  was  found  to  result  in 
better  agreement  than  the  generally  used  linear  relation 


fi  =  0.005  e 


Entrainment  Fraction: 


(4.35) 


Although  a  constant  ratio  of  mass  flux  between  the  core  and  liquid  film 
is  not  physically  correct,  it  was  found  that  the  constant  K  in  equation 
(4.3)  is  a  weak  function  of  the  film  thickness  and  flow  rates.  A  value 
of  K  =  0.95  was  used  in  equation  (4.5)  to  obtain  the  entrainment  frac¬ 
tion. 

Eddy  Diffusivity: 

The  following  non-dimensional  turbulent  eddy  diffusivity  was  used  in  the 
present  model. 


=  0.16  (y+  (1  -  JU)) 
vi  yo  \ 


1  -  e 


where 


(j)  (— )  (Wn 

D  2 


v  •  —  (V/2 
^  2 


where  is  the  single-phase  friction  factor  for  the  simulating  liquid 
alone  and  is  discussed  in  2.5.4. 


4.5  COMPARISON  WITH  EXPERIMENTAL  DATA 

The  calculated  two-phase  multipliers  and  the  experimental  data  for 
L-45/water  flow  with  annular  flow  mixing  section  are  shown  in  Figure 
4.2.  The  calculated  void  fractions  and  the  experimental  data  for  the 
same  flow  conditions  are  shown  in  Figure  4.3. 

As  expected,  the  predictions  compare  well  with  the  experimental  data 
since  the  empirical  relations  were  developed  from  the  same  data  base. 
However,  the  model  does  not  result  in  good  agreement  with  pressure  drop 
data  for  nitrogen/water  tests  of  Ref.  (52).  The  interfacial  friction 
and  entrainment  rates  are  empirically  developed  from  the  experiments 
with  two  equal  density  liquids.  As  will  be  discussed  in  Section  5,  the 
method  of  simulating  reduced  gravities  with  two  liquids  for  annular  flow 
conditions  is  questionable.  Although  the  model  provides  a  general 
approach  for  predicting  annular  flow  pressure  drop  and  void  fraction 
under  reduced  gravities,  the  empirical  relations  should  be  modified  for 
use  with  gas-liquid  flow  in  the  absence  of  gravity. 


Experimental  Values 
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Figure  4.2  Comparison  of  the  Predicted  Values  for  the  Two-Phase  Friction 
Multiplier  With  Experimental  Values. 
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Figure  4.3  Comparison  of  the  Predicted  Values  for  the  Void  fraction  With 
Experimental  Values. 


An  attempt  was  made  to  modify  the  empirical  relations.  Density  ratio 
was  introduced  in  the  interfacial  friction  relation  to  make  it  appli¬ 
cable  to  different  fluid  combinations.  The  same  power  of  density  ratio 
as  used  in  Ref.  (69) ,  equation  (4.34),  was  used.  This  modification  made 
the  predictions  more  reasonable,  but  the  entrainment  rates  were  con¬ 
siderably  higher  than  expected  and  the  resulting  two-phase  multipliers 
were  higher  than  the  measured  values  by  a  factor  of  two.  This  result 
should  have  been  expected,  since  the  correlation  for  entrainment  rate  is 
based  on  the  data  for  two  liquids  which  have  a  smaller  interfacial  ten¬ 
sion  than  gas/liquid  flows.  Previous  studies  have  shown  that  entrain¬ 
ment  rate  increases  with  decreasing  surface  tension.  In  addition, 
visual  observation  of  our  experiment  had  shown  considerably  higher 
entrainment  than  expected  with  gas/liquid  systems.  One  could  reduce  the 
entrainment  rate  by  the  ratio  of  interfacial  tensions  or  vary  the 
constant  K  in  equation  4.5  to  improve  the  agreement  with  reduced  gravity 
data.  However,  the  limited  amount  of  data  does  not  justify  such  a 
modification  at  this  point.  The  above  approach  with  the  density  correc¬ 
tion  introduced  in  the  interfacial  friction  relation  shows  promising 
trends,  but  more  reduced  gravity  gas/liquid  data  is  needed  to  develop  a 
useful  correlation  for  entrainment. 
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Section  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


In  order  to  develop  a  better  understanding  of  two-phase  flow  beha¬ 
vior  under  reduced  gravity  conditions  and  generate  a  data  base  for 
two-phase  pressure  drop  and  void-quality  relation,  an  experimental 
program  was  developed  to  simulate  these  conditions  at  earth  gra¬ 
vity.  These  tests  provide  data  for  long  duration  steady  state  fully 
developed  two-phase  flow  conditions. 

The  simulation  of  zero  gravity  by  using  two  equal  density  liquids 
works  well  if  care  is  taken  in  selecting  the  liquids  with  proper¬ 
ties  which  represent  the  gas/liquid  system.  The  wetting  capability 
of  the  liquids  as  well  as  the  relative  viscosities  and  surface  ten¬ 
sions  are  particularly  important. 

A  test  matrix  was  selected  to  result  in  turbulent  flow  conditions 
which  have  more  practical  applications  than  laminar  flow.  The  simu¬ 
lation  of  zero  gravity  two-phase  flow  under  turbulent  conditions  is 
more  appropriate  for  bubbly  flow.  This  is  due  to  the  fact  that  in  a 
two  liquid  annular  flow  condition,  the  liquid  representing  the  gas 
phase  travels  at  a  velocity  which  is  not  greatly  different  from  the 
velocity  of  the  other  phase.  However  in  a  gas/liquid  annular  flow, 
the  phases  travel  at  considerably  different  velocities  due  to  den¬ 
sity  differences.  Therefore,  the  two  liquid  annular  flow  simula¬ 
tion  can  only  provide  a  qualitative  representation  of  the  actual 
condition  by  eliminating  the  buoyancy  effects. 

The  test  results  showed  that  the  pressure  drop  and  void-quality 
relation  can  be  closely  predicted  by  the  Homogeneous  Equilibrium 
Model  (HEM).  However,  it  is  believed  that  such  a  result  is  due  to 
the  equal  densities  of  the  fluids  rather  than  an  inherent  condition 


1 

\ 
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in  the  absence  of  buoyancy.  The  two  types  of  mixing  sections  pro¬ 
duced  the  same  void  fraction  and  pressure  drops  which  means  that 
for  fully  established  flow  the  method  of  introducing  the  phases 
does  not  affect  the  flow  patterns. 

5.  The  model  developed  for  adiabatic  bubbly  flow  resulted  in  two-phase 
friction  multipliers  which  were  considerably  larger  than  earth  gra¬ 
vity  models.  There  is  no  data  presently  available  for  bubbly  flow 
at  zero  gravity.  However,  the  data  from  airplane  trajectory  tests 
for  other  flow  patterns  have  shown  approximately  80 %  to  120%  higher 
pressure  drops  than  equivalent  conditions  at  earth  gravity.  The 
model  developed  for  bubbly  flow  resulted  in  pressure  drops  which 
were  10%  to  120%  higher  than  predicted  by  HEM.  The  void  fraction 
predicted  by  this  model  is  very  close  to  HEM  prediction. 

6.  The  model  proposed  for  bubbly  flow  with  nucleation  requires  the 

critical  heat  flux  for  evaluation  of  the  pressure  drop  and  void 
fraction.  This  model  results  in  higher  two-phase  friction 
multipliers  than  predicted  under  a  similar  adiabatic  condition. 
However,  the  increase  would  only  be  significant  at  high  heat  fluxes 
which  are  unlikely  as  long  as  bubbly  flow  exists. 

7.  The  annular  flow  model  was  based  on  a  phenomenological  approach  and 

results  in  a  good  agreement  with  the  data.  However,  since  the 

empirical  relations  used  for  the  interfacial  friction  and  the 

entrainment  rate  are  based  on  the  data  from  the  equal  density 
liquids,  they  can  not  be  extended  to  gas/liquid  flow  in  zero  gra¬ 
vity.  The  model  provides  a  general  approach  for  predicting  annular 
flow  pressure  drop  and  void  fraction.  Attempts  to  generalize  the 
empirical  relations  for  gas/liquid  flow  at  reduced  gravities  showed 
promising  trends  but  more  tests  are  needed  to  develop  relations  for 
the  ':nterfacial  friction  and  entrainment  rates. 

8.  Study  of  the  convective  boiling  phenomenon  has  shown  that  the 

mechanism  of  two-phase  heat  transfer  should  not  be  affected  by  the 
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gravity  level.  For  equal  mass  flow  rate  and  void  fraction,  the 
earth  gravity  results  for  heat  transfer  coefficient  can  be  used. 
However,  since  the  available  correlations  are  generally  expressed 
in  terms  of  mass  flux  and  quality,  a  void-quality  relation  from  the 
earth  gravity  data  is  implicitly  utilized  in  these  models. 
Therefore,  if  the  heat  transfer  rates  for  the  given  mass  flux  and 
void  fraction  are  not  available,  the  correlations  developed  from 
the  earth  gravity  data  should  be  modified  to  account  for  differen¬ 
ces  in  the  void  fraction. 

Based  on  the  results  of  the  present  study  and  the  review  of  the  existing 
and  ongoing  efforts  in  reduced  gravity  two-phase  flow,  the  following 
recommendations  for  future  studies  are  made: 

o  The  bubbly  flow  model  should  be  evaluated  by  comparing  the  model  to 
data  from  zero  gravity  gas/liquid  flow.  The  empirical  parameters  in 
this  model  are  based  on  limited  data  and  should  be  modified,  if 
needed,  when  data  from  true  zero  gravity  tests  for  different  fluids 
become  available. 

o  The  empirical  relations  for  the  interfacial  friction  and  entrain¬ 
ment  rate  in  the  annular  flow  model  should  be  re-examined  when  more 
zero  gravity  data  becomes  available. 

o  There  is  a  great  need  for  experimental  data  with  adiabatic 
gas/liquid  two-phase  flow  to  establish  the  basic  parameters  needed 
for  engineering  designs.  Obviously,  fully  developed  long  duration 
tests  will  be  preferable,  but  airplane  trajectory  tests  can  also  be 
used  in  the  absence  of  such  data.  It  should  be  noted  that  a  few 
isolated  points  will  not  be  sufficient,  and  systematic  testing  for 
a  wide  range  of  mass  flow  rates  and  qualities  with  tests  for  con¬ 
firming  the  repeatability  of  the  measurements  is  needed.  This  data 
is  particularly  important  for  establishing  the  flow  regime  maps  and 
developing  the  constitutive  relations  for  annular  flow  conditions. 
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The  heat  transfer  correlations  developed  for  earth  gravity  con¬ 
ditions  should  be  modified  and  expressed  in  terms  of  void  fraction 
rather  than  flow  quality.  These  correlations  can  then  be  used  for 
reduced  gravity  conditions. 

The  boundaries  of  the  boiling  regimes,  i.e.,  the  Net  Vapor 
Generation  (NVG)  point,  Departure  from  Nucleate  Boiling  (DNB  or 
critical  heat  flux),  and  the  minimum  film  boiling  temperature  are 
dependent  on  the  gravity  level.  There  is  no  useful  data  or  models 
for  prediction  of  these  parameters  in  reduced  gravity  conditions. 
Knowledge  of  the  heat  transfer  boundaries  is  of  great  importance 
for  design  of  any  system  which  operates  under  convective  boiling 
conditions . 
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APPENDIX 


WORKING  PRINCIPLES  OF  THE  AUBURN  MODEL  1090  CAPACITANCE  PROBE 

A.  1  SYSTEM  DESCRIPTION 

The  Auburn  Model  1090  capacitance  probe  is  designed  to  measure  the  void 
fraction  for  two-phase  flow  of  fluids  with  considerably  different 
dielectric  constants.  The  probe  system  consists  of  two  major  com¬ 
ponents:  the  sensor  spool-piece  and  the  electronics  unit. 

The  sensor  spool -piece  has  an  array  of  six  electrodes  spaced  equally 
around  its  inner  circumference  and  is  installed  as  a  part  of  the  t-est 
section.  A  flow  liner,  or  protective  coating,  separates  the  electrodes 
from  the  flow  and  provides  an  unobstructed  flow  path  through  the  sensing 
volume. 

The  electronics  unit  provides  excitation  to  the  sensor  electrodes  and 
the  average  dielectric  constant,  or  capacitance  which  is  related  to  the 
void  fraction. 

A. 2  FUNCTIONAL  DESCRIPTION 

The  Auburn  1090  operates  by  continuously  measuring  the  relative 
dielectric  constant  of  the  flow  through  the  sensor.  The  dielectric 
constant  of  the  two-phase  mixture  is  a  function  of  the  dielectric 
constants  of  the  two  components  and  the  volumetric  ratio  (void 
fraction),  and  can  be  written  as 

eA  =  (1  ~  av)  +  eGav 

where 


eA  =  the  average  electric  constant 


cf  =  the  dielectric  constant  of  the  simulating  liquid  phase 

=  the  dielectric  constant  of  water  (the  simulating  gas  phase) 

«y  =  void  fraction  as  defined  in  (4.7) 

In  order  to  have  an  accurate,  uniform  measurement,  an  electrically 
rotating  sensor  field  is  provided  by  sequentially  switching  the  excita¬ 
tion  from  electrode  to  electrode.  The  opposing  three  electrodes,  the 
receive  electrodes,  are  switched  concurrently  (see  Figure  A.l).  The 
switching  rates  provide  electric  field  rotation  at  greater  than  200  Hz, 
thus  insuring  accurate  response  in  less  than  5  milliseconds. 

The  excitation  voltage  applied  on  the  sensor  electrodes  is  provided  by  a 
high  frequency  oscillator.  The  measured  current  from  receiving  elec¬ 
trodes  is  amplified  and  processed  by  the  electronics  unit.  The  basic 
operational  block  diagram  is  shown  in  Figure  A. 2. 


